Quantum quench dynamics of some exactly solvable models in one dimension 
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The dynamics of the Luttinger model and the sine-Gordon model (at the Luther-Emery point as 
well as in the semiclassical approximation) after a quantum quench is studied. We compute in detail 
one and two-point correlation functions for different types of quenches: from a non-interacting to an 
interacting Luttinger model and vice-versa, and from the gapped to the gapless phase of the sine- 
Gordon model and vice-versa. A progressive destruction of the Fermi gas features in the momentum 
distribution is found in the case of a quench into an interacting state in the Luttinger model. 
The critical exponents for spatial correlations are also found to be different from their equilibrium 
values. Correlations following a quench of the sine-Gordon model from the gapped to the gapless 
phase are found in agreement with the predictions of Calabrese and Cardy [Phys. Rev. Lett. 96 
136801 (2006)]. However, correlations following a quench from the gapped to the gapless phase at 
the Luther-Emery and the semi-classical limit exhibit a somewhat different behavior, which may 
indicate a break-down of the semiclassical approximation or a qualitative change in the behavior 
of correlations as one moves away from the Luther-Emergy point. In all cases, we find that the 
correlations at infinite times after the quench are well described by a generalized Gibbs ensemble 
[M. Rigol et al. Phys. Rev. Lett. 98, 050405 (2007)], which assigns a momentum dependent 
temperature to each eigenmode. 
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I. INTRODUCTION 

Most of the theoretical effort in the field of strongly 
correlated quantum systems over the past few decades 
has focused on understanding the equilibrium properties 
of these fascinating systems. For instance, achieving a 
complete understanding of the phase diagram of rather 
"simple" models like the two-dimensional fermionic Hub- 
bard model still remains a huge challenge. Nevertheless, 
however important these endeavors are, understanding 
the phase diagram and the equilibrium properties of the 
phases of strongly correlated systems will not certainly 
exhaust the possibilities for finding new and surprising 
phenomena in these complex systems, especially out of 
equilibrium. 

In classical systems, the existence of steady states out 
of equilibrium is well known. Very often, however, the 
properties of such states have very little to do with 
the equilibrium properties of the systems where they 
occur. Moreover, also very often their existence can- 
not be inferred from any previous knowledge about the 
equilibrium phase diagram: They are emergent phenom- 
ena. One good example is provided by the appearance 
of Rayleigh-Benard convection cells when fluid layer is 
driven out of equilibrium by the application of a temper- 
ature gradient. Although dissipation plays an important 
role in the formation of those classical non-equilibrium 
states, one may also wonder if also non-quilibrium steady 
states can appear when quantum systems are driven out 
of equilibrium. Differently from classical systems, dis- 



sipation in quantum systems causes decoherence, which 
usually destroys any quantum interference effects that 
could lead to quantum non-equilibrium states without 
classical analog. However, it is known that decoherence, 
due to coupling with the environment, is always present 
in most quantum many-body systems for large enough 
collections of particles. This may explain why the study 
of non-equilibrium phenomena in quantum many-body 
systems has been regarded, until very recently, as a sub- 
ject of mostly academic interest. 

The recent availability of highly controllable systems 
such as mesoscopic heterostructures (e.g. quan- 
tum dots) and especially ultra-cold atomic gases 
has finally provided the largely lacking experimen- 
tal motivation for the study of non-equilibrium 
phenomena, leading to an explosion of theoretical 
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Cold atomic gases are especially interesting because 
they are very weakly coupled to the environment, thus 
remaining fully quantum coherent for fairly long times 
(compared to the typical duration of an experiment). 
At the same time, it is relatively easy to measure the 
coherent evolution in time of observables such as the 
density or the momentum distribution. Thus, theorists 
can now begin to pose questions such as: Assuming that 
a many-body system is prepared in a given initial state 
that is not an eigenstate of the Hamiltonian, how will 
it evolve in time? And, more specifically, will it reach a 
stationary or quasi-stationary state? If so, what will be 
the properties of such a state? How much memory will 
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the system retain of its initial conditions? 

From another point of view, the problem described in 
the previous paragraph can be formulated as the study 
of the response of a system to a sudden perturbation in 
which the Hamiltonian is changed over a time scale much 
shorter than any other characteristic time scales of the 
system. In what follows, we shall refer to this type of ex- 
periment as a quantum quench. Quantum quenches arc 
also of particular interest to the 'quantum engineering' 
program for cold atomic gases4^ The reason is that, if 
we intend to use these highly tunable and controllable 
systems as quantum simulators of models of many-body 
Physics (such as the 2D fermionic Hubbard model men- 
tioned above), it is utterly important to understand to 
what extent the final state of the quantum simulator de- 
pends on the state in which it was initially prepared. In 
particular, one is interested in finding out whether the 
observables in the final state state can be obtained from 
a standard statistical ensemble (say, the microcanonical, 
or the canonical ensemble at an effective temperature) . If 
this is so, one would speak of thermalization. If this does 
not happen, then how much memory does the system 
retain about its initial state beyond the average energy 
E = (H)7 

We would like to emphasize that the above questions 
are not a merely academic. Indeed, cold atomic systems 
allow for the study of non adiabatic dynamics when the 
system is driven between two quantum phases such as a 
superfluid and a Mott insulator £Ldl Also, in a recent 
experiment^ it has been shown that a faithful real- 
ization of the Lieb-Liniger model 4 - 4 , exhibits absence of 
thermalization. In other words, when prepared in an 
non-equilibrium state, the experimental system reached 
a steady state that cannot be described by any of the 
'standard' ensembles of Statistical Mechanics. This ab- 
sence of thermalization seems to be a consequence of the 
integrability of the Lieb-Liniger model, that is, the exis- 
tence of an infinite number of independently conserved 
quantities. This conclusion was backed by the theoreti- 
cal analysis of Rigol and coworkers^ who observed that 
the non-equilibrium dynamics of an integrable system is 
highly constrained. Thus, based on numerical simula- 
tions for the Tonks- Girardeau limit of the Lieb-Liniger 
model, these authors conjectured that the long-time val- 
ues of some observables should converge to those ob- 
tained from a Generalized Gibbs ensemble, which can be 
constructed from a maximum entropy approach ] 10 ' 45 i 46 
The conjecture was first analytically confirmed by ana- 
lyzing an interaction quench in the Luttinger model by 
one of us^ Later, it has been also found true in other in- 
tegrable models: Cardy and Calabresse studied a quench 
in a Harmonic chain^ Eckstein and Kollar analyzed the 
Falikov-Kimball model in infinite dimensions^ and the 
1/r Hubbard model in one dimension^ Moeckel and 
Kehreir* 1 ^ studied an interaction quench in the Hubbard 
model in infinite dimensions by a flow equation method, 
and found that the system reaches an intermediate non- 
thermal state. Finally, recent numerical studies also have 



suggested that lack of thermalization may even persist in 
the absence of integrability in one-dimensional systems^ 
or that it may occur only certain parameter regimes of 
non-integrable models^ 

However, it can be expected that^ for a rather gen- 
eral choice of the initial state, along with a situation 
where there are few conserved quantities, the system will 
lose memory of most of the details of the initial state 
and, after it reaches a steady state, the expectation of 
most experimentally accessible observables such as the 
particle density or the momentum distribution, will look 
essentially as those obtained from a standard thermal 
ensemble^ Indeed, this is what seems to be observed in 
the vast majority of the experiments with cold atomic 
gases. However, for experimental many-body systems in 
general, it is hard to quantify whether this will be always 
the case. We should take into account that (except per- 
haps in the case of cold atomic gases) the exact form of 
the quantum Hamiltonian is frequently not known with 
accuracy. And even when it is known, it is not always 
possible to tell a priori whether the system is integrable 
or even if it has other conserved quantities besides the 
ones assumed by the standard thermodynamic ensem- 
bles. As a possible experimental check, we can say that, 
provided the final result is largely independent of the 
particular details of the preparation of the initial state 
and provided that some simple observables can be com- 
puted assuming that the final steady state is thermal, we 
can say that thermalization has occurred. Indeed, some 
recent numerical evidence^ supplemented by the exten- 
sion to many-particle systems of a conjecture known as 
'eigenstate thermalization hypothesis' (first introduced in 
the context of quantum chaos^i) , seems to indicate that 
lack of integrability will in general lead to thermalization 
(in the sense defined above). Indeed, Reimanr*^ has re- 
cently analytically demonstrated that, under realistic ex- 
perimental conditions, equilibration will be observed in 
an isolated system that has been initially prepared in an 
non-equilibrium mixed state. Nevertheless, even of the 
issue of thermalization for non-integrable systems may 
have been finally settled, other questions such as the de- 
tails of the transition from the integrable case (which 
thermalized to a generalized Gibbs ensemble) to the non 
integrable case (which thermalizes to the standard micro- 
canonical, or for large enough systems, Gibbs ensembles) 
are questions that are still far from being completely un- 
derstood^ 

In this article, we will not try to answer the difficult 
questions posed in the previous paragraph. Instead, we 
focus on analyzing the quench dynamics of two relatively 
well-known one dimensional models: the Luttinger model 
(for which a brief account of the results has been already 
published elsewhere^) and the sine-Gordon model. For 
the latter, we present results in two solvable limits in 
which the Hamiltonian can be reduced to a quadratic 
form of creation and destruction operators. The long 
time behavior of the single particle Green's function after 
a quench for general quadratic Hamiltonians of this form 
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has been recently studied by Barthel and Scholwock. 48 
These authors provided some general conditions for the 
appearance of dephasing and steady non-thermal states. 
This question as been also taken up recently by Kollar 
and Eckstein. 18 However, since the models taken up in 
this article may be relevant for the experiments with cold 
atomic gases (see Sect. IVI[) or numerical simulations, it 
is important to obtain analytical results for them. The 
simplicity of these models also allows us to test in detail 
a number of general results i 12 ' 48 

The rest of this article is organized as follows. In 
Sect. |TT] we describe the how the dynamics after a quan- 
tum quench can be obtained as a time-dependent canon- 
ical transformation for fairly general quadratic Hamil- 
tonians. The solution is used extensively in Sects. IIIII 
and IIVI to obtain the evolution of observables and corre- 
lation functions of the Luttinger and sine-Gordon models. 
For the Luttinger model we consider the case where the 
interaction between the fermions is suddenly switched on 
(for which some of the results where briefly reported in 
Ref. 1) and the reverse situation, when the interaction 
in suddenly switched off. For the sine-Gordon model, 
we study the situations when the system is quenched be- 
tween the gapped and the gapless phases of the model, 
and viceversa. We show that in the latter case (i.e. when 
the system is quenched from the gapless to the gapped 
phase), the results obtained at the Luther-Emery point 
(where the model can be exactly mapped to a quadratic 
fermionic Hamiltonian) and in the so-called semiclassical 
approximation are qualitatively different. The origin of 
this difference is still not well understood, and we cannot 
discard that it is indeed an artifact of the quasi-classical 
approximation (which, however, yields results in agree- 
ment with those obtained at the Luther-Emery point for 
the case when the system is quenched form the gapped to 
the gapless phase). In Sect. [V] we show the long time be- 
havior of some of the correlation functions in these mod- 
els can be obtained from a generalized Gibbs ensemble 
(along with instances where it fails). The experimental 
relevance of our results is briefly discussed in Sect. IVI1 
along with other conclusions of this work. Finally, the 
details of some of the most lengthy calculations are pro- 
vided in the appendices. 

II. QUADRATIC HAMILTONIANS 

As an illustrative calculation, let us first study the case 
of a quantum quench in a model described by a quadratic 
Hamiltonian: 

H(t) = H [(jo® + m(q, t)} b\q)b(q) 
q 

+ 1 -Y d fig{ q ,t) [b{ q )b{-q)+b\q)b\-q)] , (1) 
9 

where [b(q), tf(q')] = S q , q ', commuting otherwise. We 
will assume that the quench takes place at t = 0, so that, 



within the sudden approximation, the system is described 
by Hi = H(t < 0) for t < and by H { = H(t > 0) for 
t > 0. Furthermore, in order to simplify the analysis, 
we assume that m(q,t ^ 0) = g(q,t ^ 0) = 0, and 
m(q,t > 0) = m(q) and g(q,t) = g(q). Notice that the 
initial Hamiltonian is diagonal in the ^-operators: 

H i = H Q = Y, fiwo(g)6 f (?)6(?). (2) 

9 

In order to obtain the time evolution of operators O = 
0[{b^(q),b(q)}} after the quench, we recall that, in the 
Heisenberg picture, 0(t > 0) = e iH t t/hQ e -iH t tjH = 
0({b(q,t),b' f (q,t)}), and therefore all that is needed to 
solve the above quench problem is to obtain the time 
evolution of b(q) for t > 0. For Hamiltonians like ([T]) 
this can be done exactly because Hf = H(t > 0) can 
be diagonalized by means of the canonical Bogoliubov 
( "squeezing" ) transformation: 

a(q) = cosh/%) b(q) + sinh /?(<?) b\-q). (3) 

Upon choosing 

tanh2/3(g)= (4) 

the Hamiltonian at t > is rendered diagonal: 

Hf = H = E + J2 fiu (9)aHq)a(q), (5) 

9 

where Eq is the energy of the ground state of H (relative 
to the ground state energy of Ho) and 

u(q) = ^[^(q)+ m {q)f-[ g {q)Y (6) 

the dispersion of the excitations about the ground state 
of Hf. The evolution of the a(q) is given by a(q,t) — 
e iH{t l h a(q)e~ lHtt / h = e'^* a{q) . By application of a 
direct and reverse Bogoliubov transformation, one can 
obtain the time evolution of b(q): 

b(q,t) = f(q,t)b(q) + g*(q,t)tf{-q), (7) 

where 

f(q,t) = cosiv(q)t — is'mco(q)t cosh2/3(g), (8) 

g(q,t) = isinuj(q)t sinh2/3(q). (9) 

It is easy to check that ([7]) obeys the initial condition, 
b(q,t = 0) = b(q), and also respects the equal-time com- 
mutation rules, 

[b(q,t),b(q',t)} = (f(q,t)g*(q,t) 

-g*(q,t)f{q,t))5 q ,- q = Q, (10) 
[b{q,t),b^,t)] = (\f(q,t)\ 2 -\g(q,t)\ 2 ) 6 q , q , 

= <W (ii) 
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Thus, a quantum quench described by a quadratic Hamil- 
tonian can be solved by means of a time-dependent 
canonical transformation. 

When the quench is reversed, i.e. when the case with 
m(q,t ^ 0) = g(q,t ^ 0) = 0, and m(q,t < 0) = m(q) 
and g(q,t < 0) = g(q) is considered, the roles played by 
the initial and final Hamiltonians arc also reversed: the 
final Hamiltonian is now diagonal in the 6's, Hi = Hq, 
whereas the transformation of Eq. ([3]) renders the initial 
Hamiltonian, H = H, diagonal. Therefore, in this case 
the evolution of the 6-operators is trivial: Hf. b(q, 1) — 
e -iuj {q)t^ w h ereas the evolution of the a's is given by 

a(q,t) =Mq,t)a{q)+gZ(q,t)a)(-q), (12) 

where 

fo{q,t) = cos uJo(q)t — i sinuio(q)t cosh2/3(g), (13) 

9o{q,t) = -ismu (q)t sinh2/3(g). (14) 



III. THE LUTTINGER MODEL 

The Luttinger model (LM) is a one-dimensional (ID) 
system of interacting fermions with linear dispersion. It 
was first described by Luttinger— but its complete so- 
lution was later obtained by Mattis and Lieb^ who 
showed that the elementary excitations of the system 
are not fermionic quasi-particles. Instead, they intro- 
duced a set of bosonic operators describing collective 
density modes (phonons) of the system, which are the 
true elementary low-energy excitations of the model. The 
methods of Mattis and Lieb bear strong resemblance to 
the early work of Tomonaga 51 on the one-dimensional 
electron gas. Extending the work of Tomonaga, Mattis 
and Lieb, Luther and Peschel 52 computed the equilib- 
rium one and two-particle correlation functions, showing 
that all correlation functions exhibit (at zero tempera- 
ture) a non-universal power-law behavior at long dis- 
tances, which signals the absence of long-range order. 
Later, Haldan o 53 i 54 ' 55 conjectured that these properties 
(i.e. collective elementary excitations exhausting the 
low-energy part of the spectrum as well as power-law 
correlations) are a distinctive features of a large class 
of gapless interacting one-dimensional systems, which he 
termed (Tomonaga-) 'Luttinger liquids'. Thus, the Lut- 
tinger model can be understood as a fixed point of the 
renormalization-group for a large class of gapless many- 
body systems in one dimension. Therefore, the thermal 
equilibrium properties of many ID system are univer- 
sal in the sense that they can be accurately described 
by the Luttinger model. However, in this work we shall 
be concerned with the non-equilibrium properties of the 
LM, and because non-equilibrum phenomena can involve 
highly excited states, we shall make no claim for univer- 
sality. The precise conditions conditions under which the 
results obtained here apply to real systems that are the 
Tomonaga-Luttinger liquids should be investigated care- 



fully for each particular system (see also discussion in 
Sect. ED- 

The Hamiltonian of the Luttinger model (LM) can be 
written as follows: 

H hM = H Q + H 2 + H^ (15) 

H Q= ^2 hv FP ■ ^(P)^a(p) : , (16) 
p,a—r,l 

H2 = ^Y^g 2 (q):J r (q)J l (q):, (17) 

Hi = ^ 9i(q) ■ J a (q)J a (~q) ■ • (18) 

q,a=r,l 

The index a = r,l refers to the chirality of the fermion 
species, which can be either right (r) or left (I) moving; 
the symbol : . . . : stands for normal ordering prescription 
for fermionic operators, which is needed to remove from 
the expectation values the infinite contributions arising 
from the ground stated The latter is a Dirac sea, namely 
state where all single-particle fermion levels with p < 
are occupied for both chiralities. This defines a stable 
ground state (at least at the non-interacting level), which 
will be denoted by |0). 

A. Bosonization solution of the LM 

In this section we briefly review the solution of the 
LM. The Hamiltonian in Eqs. (fl~5|) to (fT8]) can be written 
as a quadratic Hamiltonian in terms of a set of bosonic 
operators^ First note that the density (current) opera- 
tors J a (q) = V : 4>t{P + q)^a(p) ■ obey the following 
commutation rules: 

[J a (q), J (q')} = (j^j <W,0<W, (19) 

which can be transformed into the Heisenberg algebra of 
bosonic operators by introducing: 

b(q) = -i (J^j V2 mq)M-q) - 0(-q)Mq)\ ■ (20) 

where $(q) is the step function. Note that the q — com- 
ponents (sometimes calle zero modes) require a separate 
treatment since J Q (0) = N a is the deviation (relative to 
the ground state) in the number of fermions of chiral- 
ity a = r,l. Rather than working with N r and iVj, it is 
convenient to introduce: 

N = N r + N l J = N r -N h (21) 

which, since N r and Ni are integers, must obey the fol- 
lowing selection rule (— l) N = (— 1) when the Fermi 
fields obey anti-periodic boundary conditions (ip a (x + 
L) = —ipa(x) (L is the length of the system), and 
therefore ip a (x) = L~ l / 2 Y^ p e~ a °^e lpx 4> a (p), with p = 
2(n — i)7r/L, n being an integer, and oq — > + )i^ 
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The Hamiltonian H^m can be expressed in terms of 
the bosonic operators introduced in Eq. (|2U1) : 

#o = £ hv F \q\ b\q)b{q) + (A 2 + J 2 ) , (22) 



**g»{0) {N2 j2) 



2L 



(23) 



H, = J2 %4(?)kl 6 f (?)&(?) + ^g^(A 2 + J 2 ). (24) 

Ignoring, for the moment, the zero mode part (i.e. q = 
terms, involving J and TV), the above Hamiltonian has 
the form of Eq. (fTJ), with the following identifications: 
w (q) = v F \q\, m(q,t) = g^(q)\q\, and g(q,t) = g 2 (q)\q\, 
and it can be therefore be brought into diagonal form by 
means of the canonical transformation of Eq. ([3]). Hence, 
the Hamiltonian takes the form of Eq. © with ui(q) = 
v(q)\q\, being v(q) = {[v F + g 4 {q)} 2 - [52(g)] 2 } 172 , and 
g^O. As to the zero mode contribution: 



Hzm 



~2L~ 



ft^Vj 2 

2L J ' 



(25) 



where v N = v F +g 4 (0)+g 2 (0) and vj = «f +54(0) -52(0). 
This defines the equilibrium solution of the LM. In the 
following sections we shall be concerned with the quench 
dynamics of this model. 



B. Suddenly turning-on the interactions 

Although it is possible to solve the general quench 
problem between two interacting versions of the Lut- 
tinger model, we shall focus here on the cases where 
the interactions described by H 2 and H4 are suddenly 
switched on (this section), and (next section) switched 
off. Thus, in this section, we shall assume that we have 
made the replacements 52,4(9) — * 92 4(<Z> t) = 92,i(q)@(t) 
in Eqs. (|17ll8|) . The Hamiltonian at times t > is there- 
fore the interacting LM. In other words, in the notation 
introduced in Sect. [H] H { = H + H 2 + H A = H^ M , 
whereas the initial Hamiltonian (for t < 0) is Hi = Hq. 
However, we note that, since both zero modes, J and 
N, are conserved by Hi = H and Hf = H^m, their dy- 
namics can be factored out, and we shall assume hence- 
forth that we work within the sector of the Hilbert space 
where J = N = (this sector contains the non-ineracting 
ground state, |0)). Thus, from now on, we shall omit Hzm 
in all discussions. 

As to the initial state, we shall consider that, within 
the spirit of the sudden approximation, at t — the sys- 
tem is prepared in a Boltzmann ensemble of eigenstates 
of Hi, at a temperature T: 



where Zq = Tre~ Hi / T . We shall further assume that 
the contact with the reservoir is removed at t — 0, and 
that, after the quench, the system evolves unitarily in 
isolation. 

Whereas Eq. ([7]) defines the solution to the interaction 
quench in terms of the modes that annihilate the initial 
ground state |0), the solution itself is not particularly il- 
luminating. To gain some insight into the properties of 
the system following the quench, let us compute a few ob- 
servables. Amongst them, we first turn our attention to 
the instantaneous momentum distribution, which is the 
the Fourier transform of the one-particle density matrix: 

C^(x,t) = (e mt / h 4(x)iP r (0)e- iH */ h ) , (27) 

where (• • • )o means that the expectation value is taken 
over the ensemble described by p Q (cf. Eq. Po]) ). The 
time dependence of the operators is dictated by H{, 
as described in Section [Til Notice that, since in gen- 
eral [Hf , po] ^ 0, time translation invariance is broken, 
and the above correlation function is explicitly time- 
dependent. 

The time evolution of ip a (x) can be obtained using the 
bosonization formula for the field operatori 52 i 54 ' 56 



tp a (x) 



r) a e' 



(28) 



Po = p(t = 0) 



Zo'e 



-HJT 



(26) 



being rj r , r\i two Majorana operators (also known as Klein 
factors, which in the present case reduce to two Pauli ma- 
trices) that obey {T] a ,r]/3} = 25 a p, thus ensuring the an- 
ticommutation of the left- and right-moving Fermi fields; 
we have also introduced the index s a = 1 for a = r and 
s a = — 1 for a = I. The bosonic fields 

cf)a(x) = s a ip 0a + ^-N a + $t (x) + $ a (x), (29) 

where [N a , tpop] = i&a,f)i and, in terms of Fourier modes, 

$ a (x) = lim yf— e- qa °/ 2 e ls " qx b(s a q). (30) 

The details of the calculations of C^, r (x,t) have been 
relegated to the Appendix El In this section we will 
mainly describe the results. However, a number of re- 
marks about how the calculations were performed are 
in order before proceeding any further. We first note 
that interactions in the Luttinger model are assumed 
to be long ranged i 49 ! 50 ' 54 This can be made explicit in 
the interaction couplings by writing 52,4(9) = 52,4(9^0), 
where the length scale Ro <C L is the interaction range. 
Thus, just like system size L plays the role of a cut-off 
for 'infrared' (that is, long wave-length) divergences, the 
interaction range, R plays the role of an 'ultra-violet' 
cut-off that regulates the short-distance divergences of 
the model. The results given below were derived assum- 
ing a particular form of the interaction (or regularization 
scheme) where the Bogoliubov parameter (cf. Eq. 2]) is 
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chosen such that sinh2/3(q) = 7 e I <? ^0 1 / 2 . Furthermore, 
we replaced v(q) by v = v(0). Indeed, these approxima- 
tions are fairly similar to the ones used to compute the 
time-dependent correlation functions in equilibrium, 52 
given that the expressions that we obtain for the out- 
equilibrium correlators are fairly similar as well (see Ap- 
pendix [A] for details). This regularization scheme greatly 
simplifies the calculations while not altering in a signifi- 
cant way (except perhaps for pathological cases, like the 
Coulomb interaction, where both v(q)/\q\ and sinh2/3(g) 
are singular at q = 0) the asymptotic behavior of the 
correlators for distances much larger than Rq. 

Returning to the one-particle density matrix defined 
above in Eq. (|27p. we note that it can be written as can 
be written as the product of two factors: C^ r (x,t) — 

C^(x)h r (x,t), where C$(x) is the noninteracting one- 
particle density matrix, and thus h r (x,t) accounts for 
deviations due to the interactions. Hence, this factor- 
ization allows us to write the instantaneous momentum 
distribution function as a convolution: 



/(P,*) 



k)h r (k, t), 



(31) 



where f {0) (p) = (e hVF P/ T + l)" 1 is the Fermi-Dirac dis- 
tribution, and h r (k,t) is the Fourier transform of h r (x, t) 
in x variable. 

Before presenting the results for the expression of the 
one-particle density matrix as well as the momentum dis- 
tribution at finite temperatures, it is worth considering 
the (much simpler looking) zero-temperature expression. 
We first discuss finite-size effects. For a system of size L 
we obtain: 9 



C i)r {x\L)=Cf r {x\L) 



Ro 



d{x\L) 

d(x - 2vt\L)d(x + 2vt\L) ^ /2 



d(2vt\L)d{-2vt\L) 



(32) 



where d(z\L) = L|sin(7rz/L)| is the cord function and 
C$(x\L) = i{2Lsm[ir(x + iao)/L]y 1 (a -> 0+) the 
noninteracting one-particle density matrix. Notice that 
this result is valid only asymptotically, that is, for 
d(x\L),d(x ± 2vt) ^> R . Thus, we see that C^ r (x,t\L) 
is a periodic function of time with period equal to To = 
L/2v. This is in agreement with the general expectation 
that correlations in finite-size systems exhibit time recur- 
rences because its energy spectrum is discrete. Although 
the recurrence time generally depends on the details of 
the energy spectrum, in the present version of the LM, 
the spectrum is linear uj(q) ~ v\q\ and thus, the energy 
spacing between (non-degenerate) many-body states is 
Ao — 2irhv/L. Hence, the recurrence time To ~ 27r/Ao 
follows (the extra factor of | is explained by the so-called 
light-cone effect, see further below). The periodic be- 
havior exhibited by the one-particle density matrix (127|) 
implies that, after the quench, the system will not reach 



a stationary state with time independent properties as 
t — ► +00. A similar conclusion is reached by analyz- 
ing, e.g. the finite-size version of the density correlation 
function, 9 

C Jr (x,t\L) = (e kH!t ' h J r {x)Jr{Q)e- lHtt l h ) Q 

(l+ 7 2 )/4^ 2 7 2 /8^ 2 7 2 /8^ 2 



[d{x\L)Y 



[d(x - 2vt\L)} 2 [d{x + 2vt\L)f' 

(33) 



where 



(34) 



is the density (also referred to as current) operator in real 
space. 

However, in the thermodynamic limit L — > 00, the re- 
currence time To = L/2v — > +00, and the system does 
reach a time- independent steady state. In this limit, 
d(x\L) — > \x\, and the asymptotic form of the single par- 
ticle density matrix becomes^ 



CV, r (x,t > 0) = Gi 0) {x) 
Hence, 

h r (x, t) 



Rq 


7 2 


x 2 - (2vt) 2 


7 2 /2 


X 




(2vt) 2 





Ro 




x 2 - (2vt) 2 


x 




(2vt) 2 



7 2 /2 



(35) 



(36) 



In order to understand the evolution of the momentum 
distribution, without actually having to compute it, it 
is very useful to consider the various limits of the above 
expression, Eq. (f3"5]) . First of all, for short times such 
that 2vt -C the function h r (x,t) is asymptotically 
just a time-dependent factor^ 



(37) 



which can be interpreted as a time-dependent 'Landau 
quasi-particle' renormalization constant in an effective 
time-dependent Fermi-liquid description of the system. 
In other words, as time evolves after the quench, wc 
could imagine that the quasi-particle weight at the Fermi 
level is reduced from its initial value (Z(t = 0) = 1) to 
Z(t > 0) < 1. At zero temperature, this time-dependent 
renormalization of the quasi-particle weight reflects itself 
in a reduction of the discontinuity of the momentum dis- 
tribution f(p, t) at the Fermi level (which corresponds to 
p = in our notation). Therefore, at any finite time, the 
system behaves as as if it was a Fermi liquid and there- 
fore it keeps memory of the initial state (a non- interacting 
Fermi gas). 

Yet, for t — > +00, h r (x, t) becomes a power-law^ 



lim h r (x,t) 

t — >oo 



Ro 



(38) 
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Interestingly, this time-dependent reduction of the quasi- 
particle weight after quenching into the interacting state 
has been also found in Ref. LL9| when studying a simi- 
lar quench in the Hubbard model in the limit of infinite 
dimensions. In this case, however, the discontinuity re- 
mains finite even for t — ► +00, which is different from the 
behavior of the LM, which is known to be a non-Fermi 
liquid system in equilibrium. These non-fermi-liquid fea- 
tures also persist in the quench dynamics, as we have 
found above. 



The result of Eq. is similar to the zero temper- 
ature result in equilibrium. However, the exponent of 
C^ r (x,t — + 00) is equal to 1 + 7 2 , and, even for an in- 
finitesimal interaction (i.e. 7 <C 1), it is always larger 
than the one that governs the asymptotic ground state 
(i.e. equilibrium) correlations ! 52 ' 54 7 2 = sinh 2 2/3(0) > 
7 2 q = 2sinh 2 /3(0). The reason for the larger exponent 
can be understood from two facts: i) By the variational 
theorem, the initial state (i.e. the ground state of the 
non-interacting Hamiltonian Hq), is a complicated ex- 
cited state of the system, it) Both Ho and the Hamilto- 
nian that performs the time-evolution (Hf = H^m) are 
critical (i.e. scale free, apart from the cut-off Ro), thus, 
the system is likely to remain critical. Before going fur- 
ther in the discussion of other results, we should note 
that, in the literature on Tomonaga-Luttinger liquids, it 
is customary to introduce the dimensionless parameter 
K = e" 2 ^ ), in terms of which 7 2 = (K - K~ x fl^ 
(whereas 7 oq = (K + K^ 1 ~ 2)/2). We shall use the 
parameter K in other expressions below. 

It is worth emphasizing that the particular evolution 
of the asymptotic correlations from Fermi liquid-like at 
short times to non-Fermi liquid-like at infinite time ex- 
hibited by the one-particle density matrix, is also found 
in other correlation functions. However, the idea that the 
system 'looks like' an interacting Fermi liquid at any fi- 
nite t should not be taken too far. In this regard, we 
should note that for |a;| ^> 2vt, the prefactor of the 
term oc (2irx)~ 2 of the density correlation function (cf. 
Eq (|33p h which in equilibrium is proportional to the sys- 
tem compressibility — remains equal to (minus) unity, 
which is the value that corresponds to a non-interacting 
Fermi gas (in an interacting Fermi liquid it would devi- 
ate from one). For t — > +00 the same prefactor becomes 
(1 +7 2 ) > 1, i.e., which does not directly reveal the non- 
Fermi liquid-like behavior. However, other correlation 
functions (like the one-particle density matrix discussed 
above) do. For instance, consider the following: 

C?(x,t) = ( e 2im*(».t) e -2im*(0,t)) ) (39) 

C$(x,t) = (e^Me-^t -')), (40) 

where 4>(x) — i [4>r(x) + 4>i(x)\ and 9(x) = 
\ \4> r (x) — <pi(x)}. The spatial derivatives of <fi and 
8 are related to the (total) density and current density 



Z=l 

1 ► 

(a) 



i k 








Z(t) 


M 


(b) 








z=o 



(c) 



FIG. 1: Schematic of the time evolution of the momentum 
distribution f(p, t) at zero temperature, (a) At t — 0, the mo- 
mentum distribution is that of non-interacting fermions, with 
a discontinuity at the Fermi level (p = 0) Z — 1. (b) At t > 

the discontinuity is reduced in a power-law fashion Z(t) ~ 
_ 2 

t 7 . (c) For t — > +00 the discontinuity disappears and the 
momentum distribution exhibits a power-law singularity close 
to the Fermi level p = 0, f(p, t — » +00) = 1 — const, x \p\ J . 
However, the exponent characterizing the singularity is not 
the equilibrium exponent. 



fluctuations, respectively, in the system^ Using exactly 
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the same methods as above, we find (for L — > oo): 



Cf m \x) 
C f n \x) ' 



(Ro 


\ 2 x 2 - (2vt) 2 


\2vt 


J X 2 


'Ro\ 


2 x 2 ~ (2vt) 2 


K 2vt ) 


X 2 



, (41) 



n 2 (K- 2 -l)/S 



(42) 



where C^ m ^(x) = A^\R /x\ 2m2 and C {0 - n \x) = 
A e n \R /x\ n2 l 2 are the non-interacting correlation func- 
tion (where Af^ and A 9 n are non-universal prefactors). 
We note that the usual duality relation where <f> —> 8 
and K — > K , which one encounters when studying 
equilibrium correlation functions^ still holds for these 
non-equilibrium correlators. Let us next analyze their 
asymptotic properties. We consider only C™(x,t), as 
identical conclusions also apply to Cg(x,t) by virtue of 
the duality relations. For \x\ ^> 2vt, we have: 



C?(x,t) = C (0 ' m) (x) 



Ro 
2M 



m 2 (K 2 -l) 



(43) 



Thus, up to the time-dependent pre-factor, correlations 
take the form of a non- interacting system of Fermions, 
C(°< m )(x). However, in the opposite limit (|x| <C 2vt), 
this correlator exhibits a non-trivial power-law: 



cr(x,t) 



Ro 



2 (K 2 + 1) 



(44) 



Notice that this expression also describes the infinite- 
time behavior, which is controlled by an exponent equal 
to m 2 (K 2 + 1), being again different from the exponent 
exhibited by the same correlator in equilibrium, which 
equals 2m 2 [cosh 2/3(0) - sinh 2/3(0)] = 2m 2 K. 

In order to understand why the behavior found in the 
correlations for t — ► oo in Eqs. (|38l44p , also holds for 
|x| <C 2vt, let us consider the initial state at zero tem- 
perature, po = 1 0) (0 1 1 12 i 15 As mentioned above, this is a 
rather complicated excited state of the Hamiltonian that 
performs the time-evolution, Hi = Hlm ■ This means 
that, initially, there are a large number of excitations of 
Hi, namely, phonons with dispersion to(q) — v(q)\q\. The 
distribution of the phonons (b^ (q)b(q))o = sinh 2 /3(g) is 
time-independent and peaked at q = 0. Thus, within 
the approximation where v(q) ~ v(0) = v, the excita- 
tions propagate between two given points with velocity 
v. Thus, if we consider the correlations of two points A 
and B separated a distance |x|, the nature of the corre- 
lation at a give time t depends on wether the excitations 
found initially at, say, point A, have been able to reach 
point B or not. This is not the case if |x| > 2vt, and thus 
correlations retain essentially the properties they had in 
the initial state. Thus, up to a time-dependent prefac- 
tor, C^ r (x, t) oc ci (x). However, if the two points have 
been able to 'talk to each other' through the excitations 



present in the initial state, then correlations will be qual- 
itatively different. This happens for t = to, when the 
phonons propagating from A meet the phonons traveling 
from point B, that is for x — vto = vto, or to — x/2v 
(we assume x > without loss of generality). Thus, for 
given separation x and time t, there is a length scale 2vt, 
which marks the transition between two different regimes 
in the correlations. In the instantaneous momentum dis- 
tribution, this reflects itself in a crossover as a function 
of time from a momentum distribution n(p) exhibiting a 
discontinuous Fermi liquid-like behavior, which is valid 
i.e. for |p| < (2vt)- 1 ) to a power-law behavior of the 
form ~ |p_R | 7 "sgn(p), which applies for |p| 3> (2vt)~ 1 
but \p\ <C Ro 1 (for |p| » Ro 1 we recover the free parti- 
cle behavior corresponding to the Fermi-Dirac distribu- 
tion function at T = 0). In the t — > oo limit, by using 
the regularization scheme described above, the asymp- 
totic momentum distribution at zero temperature can be 
obtained with the help of tables^ The resulting formula 
behaves as the non-interacting Fermi-Dirac distribution 
for \p\ 3> ^o^ 1 , whereas for \p\ ^> R^ 1 it describes a non- 
Fermi liquid-like steady state: 



f(p,t -> +oo) 



1 pRo 

2 2 

+ K„2_ 3 (\pRo\)C- l 2_ 1 (\pRo\) 



K^AlpRolK^-AlpRol) 
(45) 



where K v (z), C v {z) are the modified Bessel and Struve 
functions^ respectively. This expression yields a power 
law for \pRo\ <C 1, where n(p, t — > +oo) ~ | — const, x 

(|pi?o|) 7 sign(p). Note that the momentum distribution 
n(p = 0,t) = 5) which is given by the invariance of the 
LM under particle-hole symmetry ip a (p) —> Wai - P)- 

Let us finally present the generalization of the above re- 
sults for the one-particle density matrix to finite temper- 
atures, T > 0. For T <C hv F R (but T > A = 2nhv/L, 
so that we can neglect finite-size effects and effective take 
the thermodynamic limit) C^ r (x,t) takes the following 
asymptotic form: 



C^{x,t>Q\T) = cf(x\T) 



ttRo/X 



dh(x\T) 

dh(x - 2vt\T)dh(x + 2vt\T) ^ /2 



dh(2vt\T)dh(-2vt\T) 



(46) 



where C^(x\T) and dh(x\T) can be obtained from 

C^p (x\L) and d(x\L) by replacing Ls\tl{-kx/L)/it by 
A sinh(7rx/A), where A = hvp/T is the thermal corre- 
lation length. At long times, h r (x,t\T) reduces to 



h r (x) 



ttRo/X 



sinh (ttx/X) 



(47) 



Therefore we again find that C^, r {x, t — > oo|T) has a form 
similar to the the equilibrium correlation function at fi- 
nite temperature with a different exponent controlling 
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the asymptotic exponential decay of correlations. No- 
tice that the exponential decay the correlations for t > 
is a direct consequence of the fact that the initial state 
has a characteristic correlation length, the thermal cor- 
relation length The exponential decay of correla- 
tions at finite T implies that the the steady state will 
be reached exponentially rapidly approached in a time 
of the order of h/T. It is also worth noting, however, 
that the above expression depends parametrically on the 
thermal corelation length A, and it is the Fermi veloc- 
ity, Vp, which enters in the expression for A = hvp/T, 
instead of the (renormalized) phonon velocity which en- 
ters in the thermal length A cq = hv/T, characterizing 
the equilibrium correlations. Thus, since the velocity ap- 
pears only through the definition of the thermal corre- 
lation length A, or, in other words, in combination with 
the temperature, the change from v to vp can be also 
understood as an change in the temperature scale. Fur- 
thermore, in a system with Galilean symmetry,— we have 
that vK — vp and thus the parameter that controls the 
temperature scale now is the Luttinger parameter K , so 
that the asymptotic correlations at t — > oo can be re- 
garded as the equilibrium correlations with a different ex- 
ponent and an effective temperature, T e ff = T/K. Thus, 
for repulsive interactions (i.e. K < 1) we could say that, 
besides modifying the exponent, quenching the system 
into the interacting system increases the effective temper- 
ature, whereas for attractive interactions (i.e. K > 1) the 
effective temperature is reduced after the quench. This 
effect has an impact on the momentum distribution at 
finite temperatures. To demonstrate it, we need to com- 
pute the Fourier transform of h r (x). This can be done 
by relating it to an integral representation of the asso- 
ciated Legendre function P"(z)r& and thus the Fourier 
transform of h r (x) can be written as: 




Hence, the momentum distribution can be obtained by 
numerically evaluating the convolution with the Fermi- 
Dirac distribution function (cf. Eq. IpTl]) ) of the above ex- 
pression, Eq. (j4"8|) . In Figs. f2]and[2]the momentum distri- 
bution of the interacting system in the infinite-time limit 
is displayed for a non-interacting LM that is quenched 
into an interacting state with repulsive (corresponding 
to K — 0.6) and attractive (corresponding to K = 1.7) 
interactions, respectively. 




FIG. 2: Infinite time limit of the momentum distribution vs. 
Xp (A = hvp/T is the thermal correlation length in the initial 
state) for a non-interacing Luttinger model at finite temper- 
ature T that is quenched into an interacting state with re- 
pulsive interactions (corresonding to a Luttinger parameter 
K = 0.6). 




FIG. 3: Infinite time limit of the momentum distribution vs. 
Xp (X = hvF /T is the thermal correlation length in the initial 
state) for a non-interacing Luttinger model at finite temper- 
ature T that is quenched into an interacting state with at- 
tractive interactions (corresonding to a Luttinger parameter 
if =1.7). 



C. Suddenly turning-off the interactions 

Next we briefly consider the opposite situation to the 
one considered above, namely the case where in the initial 
state the fermions interact and this interaction suddenly 
disappears. The fact that the initial state is a highly 
complicated state of the Hamiltonian that performs the 
time evolution (in this case iff = Hq, cf. Eg. [15)) . implies 
that we cannot expect that a Fermi liquid will emerge 
asymptotically at long times after the quench. Indeed, 
at zero temperature a thermodynamically large system 
approaches a steady state exhibiting equal-time corre- 
lations that decay algebraically in space. However, the 
exponents differ again from the (non-interacting) equi- 
librium ones. This can be illustrated by, e.g. computing 
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the following correlation functions: 

C™(x t) = (^ e 2imS ( x ' t ) e ~ 2ime (°^)^ 

=to (AY 



\2v F t ) 
x 2 - (2v F t) 2 



\K- Y -K)I1 



(49) 



iJo_V 



2v F t J 



x 2 - (2v F tf 



where 



IS (x) 



R 



x 



x 



2m'K 



2 /2K 



(50) 

(51) 
(52) 



are the correlation functions in the initial (interacting) 
ground state (Ag/^ are non- universal prefactors). We 
note again that the duality 9 — ► <f> and K — > K^ 1 also 
holds in this case. The correlations in the stationary state 
that is asymptotically approached at long times read: 



lim CF(x,t) = A e 

t — > + CO 



lim C${x,t) = A 4 

t—>-\-00 ^ 



R 



x 

Ro 



X 



(53) 
(54) 



However, at short times, t <C \x\/2v F , correlations look 
like those of in the initial state, up to a time-dependent 
prefactor: 



Cf\x,t < \x\/2v F ) = 



2v F t 



m 2 (K- 1 -K) 



CJ(«.*<N/2«f)=(^) 



(56) 



In this case the time-dependent prefactor has also a 
power-law form. 



IV. THE SINE-GORDON MODEL 

Let us next turn our attention towards another class 
of models in one dimension whose spectrum is not nec- 
essarily always critical (i.e. gapless) like the case of the 
Luttinger model (LM) analyzed previously. This is the 



case of the sine-Gordon model, which is described by the 
following Hamiltonian: 



H sG (t) = H - hvg<y P [ dx cos2<£, 
™o J 

Ho = ^ J dx : K- 1 (d x ^ f + K 



(57) 

(d x 6) 2 : , (58) 



where : . . . : stands for normal order of the operators, 
ao is a short-distance cut-off, and the phase and density 
fields 9(x) and <fi(x), which have been introduced in the 
connection with the LM studied in the previous section, 
are canonically conjugated in the sense that the obey: 
[4>(x),d x '6(x')] = iirS(x ~ x'). This model can be re- 
garded as a perturbation to the LM, which still yields an 
integrable model. In equilibrium, the model is known to 
have two phases, which, according to the renormalization 
group analysis 56 ' 58 and for infinitesimal and positive val- 
ues of the coupling in front of the cosine term, roughly 
correspond to K < 2 (gapped phase) and K > 2 (gapless 
phase) . 

In order to study the non-equilibrium (quench) dynam- 
ics, we will consider two different types of quenches: the 
quench from the gapless to the gapped phase and the 
reversed process, from the gapped to the gapless . In 
the first case, we assume that the dimensionless coupling 
g(t) is suddenly turned on, i.e. g(t) — g6(—t). With 
this choice, = H S Q(t < 0) is a Hamiltonian whose 
ground state exhibits a frequency gap, m, to all excita- 
tions, whereas H{ = H s Q(t > 0) has gapless excitations. 
Conversely, in the second case, we consider that g(t) is 
suddenly turned off, i.e. g(t) — g0(t). In this case, the 
ground state of Hi is gapless whereas the Hamiltonian 
performing the time evolution, iff, has gapped excita- 
tions. However, although both Hi and Hf define inte- 
grable field theories, for a general choice of the parame- 
ters K and g, the quench dynamics cannot be analyzed, 
in general, by the methods discussed above. Neverthe- 
less, in two limits, the Luther-Emery point (which corre- 
sponds to K = 1, see Sect. UVXf and in the semiclassical 
limit (that is, for K -C 1, Sect. IIVBI) . it is possible to 
study the quench dynamics by the methods of Sect. HI1 
However, the statistics of the elementary excitations hap- 
pens to be different in these two cases. 



A. The Luther-Emery point 



Let us start by considering the sine-Gordon model, 
Eq. (|59p . for K = 1, which is the so-called Luther-Emery 
point. It is convenient to introduce rescaled density and 
phase fields, which will be denoted as ip(x) — K~ x l 2 (\>{x) 
and (p(x) = K 1 / 2 9(x). Thus, the Hamiltonian in Eq. 
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(l57l) becomes: 



H, 



sG 



cfc cos K(p, (59) 



where k = 2\f~K. At the Luther-Emery point k = 2 
(i.e. if = 1) and the model can be rewritten as a one- 
dimensional model of massive Dirac fermions with mass 
by using the bosonization formula for the Fermi field op- 
erators, Eq. (j2"5j) . To this end, we set <j) r {x) = tp{x) + <p{x) 
and (j>i(x) = <p(x) — <f(x). Furthermore, for computa- 
tional convenience, we choose the Major ana fermions in 
Eq. (|28p to be r\ r = <j x and r\\ — ia y . In addition, we 
note that the gradient terms in Eq. (|59p can be written 
as the kinetic energy of free massless Dirac fermions in 
one dimension i 54 i 56 i 58 

H = -ihv / dx : ipl(x)d x ip r (x) - il>}(x)d x i(>i{x) : (60) 



As far as the cosine operator is concerned, the bosoniza- 
tion formula, Eq. (|28|). implies that: 



(61) 
(62) 



if)l(x)i/}i(x) + ipUx)if! r (x) = cos2(p(x) 

ira 

r 

= — : cos2(p(x) : 



where T = ia x a y . This is almost the cosine term 
in the sine-Gordon model (cf. Eq. 159"]) except for the 
presence of the operator T. However, we note that 
r 2 = 1 and that this operator also commutes with Hq 
and with operator in the left hand-side of Eq. (fBTj) . 
The first property implies that the eigenvalues of T are 
±1 whereas the second property implies that 7?le = 

Ho+hvg(t) J dx i/jI(x)i/ji(x) + ij}J(x)il) r (x) andTcanbe 

diagonalized simultaneously . Upon choosing the eigen- 
space where T = — 1, we obtain that 



Ale(*) = —itiM J dx : tpl(x)d x tp r (x) - ij)\(x)d x ipi(x) : 
+ hvg(t) [ Ul(x)Mx) + ^i(x)Mx)\ , (63) 



is equivalent to Eq. (fB"9"|) when k = 2. 

To gain some insight into the phases described by the 
sG model, let us first consider the Luther-Emery Hamil- 
tonian, J/le in t wo (time-independent) situations: i) 
g(t) = (the gapless free fermion phase, which coin- 
cides with the LM for K — 1), and ii) g[t) = g > 0, i.e. 
a time-independent constant (which corresponds to the 
gapped phase). In order to diagonalize the Hamiltonian 
in these cases, it is convenient to work in Fourier space 
and write the fermion field operator as: 



ip a {x) 



-ao\p\ e ipx 



(64) 



where a = r,l. The limit where the cut-off cto — > + 
should be formally taken at the end of the calculations, 
but in some cases we shall not do it in order to regularize 
certain short-distance divergences of the model. It is also 
useful to introduce a spinor whose components are the 
right and left moving fields, and which will make the 
notation more compact: 



*(p) 



1p r (p) 



«(?) = 



fuv (p) 








(65) 



Thus the Hamiltonian for the gapless phase, H , reads: 
H = J2-^(P)-no(p)-*(p):, (66) 



where loq(p) = vp is the fermion dispersion. However, 
the Hamiltonian of the gapped phase, corresponding to 
git) = g > 0, H, is not diagonal in terms of the right and 
left moving Fermi fields. In the compact spinor notation 
it reads: 



ff = 2:*t(p). W (p).*(p) 



where 



Hip) 



hm 



tvm 
-hujoip) 



(67) 



(68) 



being m = vg. Nevertheless, H can be rendered diagonal 
by means of the following unitary transformation: 



*(p) 
being 







_1pvip)_ 





cos#(p) sin#(p) 
— sin6(p) cos^(p) 



lj) r ip) 
Iplip) 



tan 26 (p) 



wo(p)' 



, (69) 



(70) 



Thus the Hamiltonian of the gapped phase, in diago- 
nal form, reads (we drop an unimportant constant that 
amounts to the ground state energy): 



h = J2 Mp) [■■ ^MMp) - i>l(p)Mp) 



(71) 



where u)(p) = ^uj^ip) 2 + to 2 . We associate tj)' v ip) (i>l(p)) 
with the creation operator for particles in the valence 
(conduction) band. 

Before considering quantum quenches, let us briefly 
discuss some of the the properties of the ground states 
of the Hamiltonians Ho and H . In what follows, these 
states will be denoted as |$o) and 1$), respectively. As 
mentioned above, the spectrum of Hq is gapless, and the 
fermion occupancies in the ground state |$o) are: 

n r (p) = (®o\4(p)Mp)\$o) = 0(-p), (72) 
mip) = ($o|V>/(p)^(p)l$o) - 6ip). (73) 
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That is, all single-particle levels with negative momen- 
tum are filled out, as described in Sect. [HT] (recall that H 
is just the Luttinger model with K = 1, i.e. gi — g± = 0). 
However, H has a gapped spectrum and, therefore, when 
constructing the ground state, |$), only the levels in 
the valence band (which have negative energy) are filled, 
whereas the levels in the conduction band remain empty: 



n v (p) = ($|V4(p)Vv(p)|$) = 1, 
n c (p) = ($|^0#c(p)|*> = 0. 

}uench from the gapped to the gapless phase 



(74) 
(75) 



The first situation we shall consider is when g(t) = 
g 0(—t) in Eq. (|63|) . so that the spectrum of the Hamil- 
tonian abruptly changes from gapped to gapless (i.e. 
quantum critical). As we did in Sect. IIII1 we denote 
Hi = H LE (t < 0) = H and H { = H LE (t > 0) = H . 
Although the expressions presented below can be com- 
puted for finite temperature, T > 0, where the initial 
state corresponds to p; = e~ Hi / T /Z\, we shall restrict 
ourselves to the zero temperature case, where the initial 
state pi = |$)(<I>|. Notice in this state, (<&| cos2^(a;)|$) = 
($|cos2<p(x)|$) =Re($|e- 2i ^ iB )|$) = -(^(x)^z(a;)) ^ 
(the minus sign stems from the eigenvalue of the op- 
erator r = r) r r]i), whereas in the ground state of Hq the 
expectation value of the same operator vanishes. There- 
fore, it behaves like an order parameter in equilibrium, 
and we can expect that it exhibits interesting dynamics 
out of equilibrium. Indeed, 



-2i<p(x,t) 



= -tZ)<^(p»*M(p.*)> 



(76) 



sin20(p), (77) 



where, in the last expression, we have already taken the 
T — > limit and set (tpj. {p)ipi (p)) = — \ sin 26(p) , as it fol- 
lows from Eqs. (|69I74I75[) . The above expression can be 
readily evaluated by recalling that sin 29 (p) = m/u)(p), 
which yields, in the L — > oo limit, 



D — 2i(p(x,t)\ 



dp cos2u>Q(p)t 



(e r ^ ) = —tu 

tn \ T , ,„ . 1 / m 



(78) 



where K$ is the modified Bessel function. Thus we see 
that the 'order parameter' (cos 2<p(0, t)) decays exponen- 
tially at long times at T = 0. The decay rate is propor- 
tional to the gap between the ground state (the initial 
state) and the first excited state of the initial Hamil- 
tonian Hi = H. The existence of this gap means, in 
particular, that correlations in the initial state between 
degrees of freedom of the system are exponentially sup- 
pressed beyond a distance of the order £ c « v/m. Since 



the system is quenched into a situation where the Hamil- 
tonian performing the time evolution is critical, that is, 
characterized by excitations that propagate along light 
cones x ± vt, the light-cone argumen t 12 i 15 discussed in 
previous sections applies and translates the correlation 
length scale in the initial state into an exponential decay 
in time of the order parameter. The exponential decay 
found in the present case (a quench from a gapped to 
a gapless or critical system) is also found in the semi- 
classical approximation to the sine-Gordon model (see 
Sect. IIV Bl below) . and it is in agreement with the results 
of Calabrese and Card y 12 i 15 obtained using a mapping to 
boundary conformal field theory (BCFT). 

Next we consider the (equal-time) two-point correla- 
tion function of the same object: 



G(x, t) = ( e -« V (*.*) e a.V(o,*)) = _L ' 

Pl,P2,P3,P4 
V „ [^0 (Pi ) + ^0 (P2 ) - 1^0 (P3 ) — WO (j>i)]t 



J(pi-P2)x 



x witoiWifaWlfa^fa)). (so) 

Applying Wick's theorem, there are three different con- 
tractions of the above four fermion expectation value, 
which can be evaluated using Eqs. (|69l74l75[) .This yields: 



(81) 



(^( P )Mp)) = ($0#l(p)) = -^sin2%), (82) 

(4(P)MP)) = sin 2 9(p), (83) 

(Mp)^i(p)) =cos 2 %). (84) 
Hence, for x ^ 0, we obtain 



/ TO V 



[K [2mt)X 



Ki 



(85) 



Let us examine the behavior of this correlation function 
in the asymptotic limit where |x| ^> ^ c ~ v/m and 2vt 3> 
Since the Bessel functions decay exponentially for 
large values of their arguments, the asymptotic behavior 
depends on whether t < \x\/2v or t > \x\/2v: 



5(x,t) 



IGttvx 

m, 

32nv 2 t t 



-2m\x\/v t>\x\/2v, 
t < \x\/2v. 



-4mt 



(86) 



These results are in also agreement with those obtained 
by Calabrese and Cardy for quantum quenches from a 
non-critical into a critical state . 12 ' 15 



2. Quench from the gapless to gapped phase 

We next consider the reversed situation to the one dis- 
cussed in the previous subsection. In this case, we set 
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g(t) = g 0(t), i.e. the initial state is critical and corre- 
sponds to the ground state of H[ — H , whereas the time 
evolution is performed according to Hf = H. To compute 
the same correlations functions as above, it is convenient 
in this case to obtain the time evolution of the opera- 
tors ip r (p) and ipi (p) , whose action on the initial state is 
known [cf. e.g. Eqs. (|72l73p ]. Once again, we restrict 
ourselves to the T = case. However, obtaining finite 
temperature results from the expressions below should 
be rather straightforward. 

We first note that the time-evolved Fermi operators 
can be related to the operators at t = by means of the 
following (time-dependent) canonical transformation: 

MP,*) = e lHit Mp)e- lHtt/h (87) 

= f(p,t)Mp)+9*(P,t)Mp), (88) 

ifnfat) = e iHtt iln(p)e- iH < t ' R (89) 

= g*(p,t)A(p)+f*(p,t)ipi(p), (90) 

where f(p,i) — cosuj(p)t — i cos 29 (p) sm to (p)t and 
g(p, t) = i sin 29 (p) smuj(p)t. It can be easily shown that 
the above transformation respects the canonical anti- 
commutation relations characteristic of Fermi statistics. 
Using Eqs. (|88|) to (f90|) . we can now compute the decay 
of the order parameter. The calculation yields: 



-2i<p{ 



x,t) > = - I E Re [/*(P'*)ff(P' < )] 



(91) 



p>0 



In deriving the above expression we have used that 
f(-P,t) = f*(p,t) and g(-p,t) = g{p,t), which fol- 
lows from cos2#(— p) = — cos29(p) because cos29(p) = 
wo(p)/o;(p). Upon setting Re[f* (p,t)g(p,t)} = 
— cos29(p) sin 29{p) sin 2 ui(p)t and taking L — > +oo, we 
find: 



dp mu>o (p) 

^Krif' 



siW 2 iv(p)t (92) 



~ A(ma ) + sm.2mt + 0(r 2 ). (93) 

ivt 

The first term is a non-universal constant that depends 
on the energy cut-off clq introduced above (cf. Eq. 
Thus, we conclude that, when quenched from the critical 
(gapless) phase into the gapped phase, the order parame- 
ter exhibits an oscillatory decay towards a non-universal 
constant value, A(mao). 

Using similar methods the (equal-time) two-point cor- 
relation function can also be obtained, yielding the fol- 
lowing result: 

Q(x,t) = te 2U P(x,t) e -2i<p(0,t)\ _ / e 2i(p(s,t)w e -2iip(0,tK 



p>0 



(94) 



where (for L — > oo) the function TL{x,t) is defined as 



/*°° dp 

H(x, t) = i — sin pa; [1 — cos 2uj(p)t] 
Jo 2lT 



Hp)] 2 



(95) 



Next we obtain the behavior of this function for t — > 
+00, in which case, the term cos2uj(p)t oscillates very 
rapidly and therefore can be dropped. Upon performing 
the momentum integral, we obtain the following result 
for large \x\: 



7i(x, t — > oo) 



2ttx 2-Kra 1 x i 



(96) 



The first term in the above expression exactly cancels the 
second in the last expression in the right hand-side of Eq. 

Thus, the two-point correlation functio in the limit 
t — * oo exhibits the following asymptotic behavior: 



lim Q(x,t) = [A(mao) 

t — >+oc 



4v 4 



(27T) 2 TO 4 X 6 ' 



(97) 



This result is clearly different from the equilibrium be- 
havior of the same correlation function in the gapped 
phase, where it decays exponentially to a constant ] 56 ! 58 
Indeed, as we have just obtained, both the order param- 
eter and its two-point correlations (Eqs. l93l and |97|) . ex- 
hibit instead an algebraic decay to constant values. 



B. The semiclassical approximation 

A controlled approximation to the sine-Gordon model 
(cf. Eq. 159")) can be obtained in the limit where k 
1 (which corresponds to small K limit in the original 
notation of Eq. I57[) . In this limit, we can expand the 
cosine term in (|59p about one of its minima, e.g. tp = 
0. Retaining only the leading quadratic term yields the 
following quadratic Hamiltonian for the boson field f{x): 



sG 



: {d x ip{*)) 2 + K (d x tp(x)f 



hvg{t)n 2 

~Z 2-k 2 /2 



dx : tp 2 {x) : (98) 



Within this approximation, the problem of studying a 
quantum quench in the sine-Gordon model becomes akin 
to the general problem studied in Sect. [TTJ To see this, 
we first introduce the expansion in Fourier components 
of ip{x), 

(p{x) = 7^ +l 7Wi 5N 

where tOo(q) — v\q\; the 6-operators introduced above 
obey the Heisenberg algebra: 



(100) 



commuting otherwise. The first two terms in Eq. (|99p are 
the so-called zero-modes, whose dynamics is only impor- 
tant at finite L. In what follows we restrict our attention 
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to the thermodynamic limit (L — > oo) and therefore ne- 
glect the dynamics of those zero modes. Introducing (j9"5|) 
into (|98|) , the Hamiltonian takes the general form of Eq. 

flJJ with g(q,t) = m(q,t) = 2vg(t)n 2 /\q\a~ R ' '' /2 . Follow- 
ing the procedure described in Sec. [Til this Hamiltonian 
can be diagonalized by means of the canonical transfor- 
mation of Eq. fl3J. Introducing m 2 = Agv 2 k 2 / 'a^' K 
which plays the role of the gap in the frequency spectrum 
of the (gapped) Hamiltonian, and setting m(q) = g(q) = 
m 2 /2uj Q (q) in Eq. Q yields: 



tanh2/3(g) 



m 2 /2 



u; 2 (q)+m 2 /2- 



(101) 



The dispersion of the excitations in the gapped phase is: 



U)(q) = y/uj (q) 2 + m 2 



(102) 



As we did when analyzing the sine-Gordon model at the 
Luther-Emery point, in what follows, we shall consider 
the cases of a quench from the gapped to gapless phase 
of the above quadratic Hamiltonian, 69 Eq. ([55)1 . and the 
reverse situation, from gapless to gapped. 



1. Quench from the gapped to the gapless phase 

Let us begin the situation where the initial state is 
the ground state of Hi = H sc and at t = the Hamil- 
tonian is changed to Hf = Hq. In this case, the evo- 
lution of the expectation value of the order parame- 
ter operator e~ 2l ^ x ^ — e~ lKip ^ or the associated cor- 
relation functions can be obtained from the knowledge 
of the two-point (equal time) correlation function out 
of equilibrium for the boson field ip(x), i.e. C(x,t) — 
(ip(x, t)(p(0, t)) — (ip 2 (0, t)). To compute this object when 
Hamiltonian changes from Hi = H (gapped spectrum) 
to Hf (gapless spectrum), we first insert the Fourier ex- 
pansion of <p(x), Eq. ([nU), and use that the time evo- 
lution (as dictated by Hi = Hq) of the ^-operators is 
b(q,t) — e~ llA)a ^ t b(q). Finally, we use the Bogoliubov 
transformation, Eq. ([3]), to relate the &-to the a-bosons, 
the basis where Hi is diagonal, and obtain the following 
expression: 



(^,tMo,t)) = -^ 



2ttv 



sinh2/3(g) 



x cos (qx - 2cu (q)t) (2n B (q) + 1) 

— 2 cosh 2p(q) cos qx Ub (q) 

- e lqx sinh 2 f3{q) - e~ lqx cosh 2 fi{q) 

(103) 

where «b(<?) = (e huJ ^ q ^ T — 1) _1 is the Bose factor (we 
have assumed that the initial state is given by the density 



matrix pi 



-Hi/T 



/Z). 



Next let us consider the behavior of the expecta- 
tion value of the order parameter. Taking into account 



that ( e -2#(o,t)) = ( e -*«vCo.*)) = e -T-<v 2 (o,*)) j we see 
that (Lp 2 (0,t)} must be evaluated in closed form us- 
ing Eq. (|103[) . Before performing any manipulation of 
this expression, it is convenient to subtract the constant 
(<p 2 (0,0)), which is formally infinite (i.e. it depends on 
the short distance cut-off, do)- Thus, (at T = and 
L — ► oo) we have: 



<^(0,t))-(^ 2 (0,0)) 



+ oo 



d(qv) 



sin 2 LO (q)t. (104) 



Inserting the expressions for oj(p) and u>o(p) in the above 
equation, we obtain: 



(^ 2 (0,i))-(^ 2 (0,0)) = -/(2mi/ft), 
where f(z) is defined as: 



3/2 

1 1/2/' 



f(z) = l + -Gfl - 



(105) 



(106) 



being G\ \ the Meijer G function^ Using the asymptotic 

expansion for this function, f(z) ml— ^Jp, and hence 
the long-time behavior of the order parameter is: 



We next examine the behavior of the two-point corre- 
lation function of the same operator, 

G(x,t) = ( e «*(*.*) e -a^(o,t)j = e i #c( a: ,t) ) ( 1Q8 ) 

where we have defined C(x,t) — (cp(x,t)cp(0,t)) — 
(ip 2 (0,t)). At zero temperature, with the help of Eq. 
(flOg)) , we find that 



C(x,t)-C(x,0) = - — J 



d(vq) 



o M<7)]M<7) 
x (1 — cosgx)(l — cos2uj (q)t), (109) 



where 



„. . 1 f°° d(vq) , 
C(x,0) = -- — i 1 ~ cosqx) i 
2 Jo u{q) 



(110) 



being clq is the short-distance cut-off. Evaluating the 
integrals: 

C(x, t) - C(x, 0) = f(mx/v) + /(2m*) 

f[m(x/v + 2t)] + f[m(x/v - 2t)] 

2 ' [ ' 



where f(z) has been defined in Eq. (|106[) . Thus, asymp- 
totically (for max{x/2w,t} ^> m _1 ): 



g(x,t) = e« <K*,0)x{__^ 



e K ^|x|/2^ )for t > \ x y 2v 

e" 

(112) 
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where Q{x, 0) describes the correlations in the initial 
(gapped ground) state, and exhibits the following asymp- 
totic behavior: 



we arrive at the following result (at zero temperature, 
and for L — > +00): 



m 2 f c 

G(x, 0) ~ B(ao) ( 1 - n 2 — e ~™MM ( U3 ) C(x, t) - C(x, 0) = — / 

\ rax J Jo 



d{vq) 



r(l — cos qx) 



where B(ag) is a non-universal prefactor. Thus we see 
that the asymptotic form of G(x, t) (Eq. lll2[) , as well that 
of the order parameter, Eq. (|107p . have the same form 
as the results found at the Luther-Emery point, and also 
agree with the results of Calabrese and Cardy based on 
BCFT^ 2 .^ 



2. Quench from the gapless to the gapped phase 

In this case, the system finds itself initally in the 
ground state of Hi = H Q , and suddenly (at t = 0) the 
Hamiltonian is changed to Hf = H sc . For this situa- 
tion convenient to obtain the evolution of the observ- 
ables from the time-dependent canonical transformation 
of Eq. ([7]), where /3(g) and tu(q) are given by Eq. QlOip 
and Eq. (|102[> . respectively. In this case, 



1 



2ttv 
Lo{q)L 



sinh 2(3(q) 



1) 



x cos (qx — 2u(q)t) (2rts(g) 
+ 2cosh2/3(q) cos qxns(q) 

+ e tqx sinh 2 f3(q) + e~ lqx cosh 2 /3(q) 

(114) 

where /ib(<?) = (e^ ^^ — 1) _1 is the Bose factor in the 
initial state (i.e. the gapless phase). 

As in the previous case, (e -2 ^^)) = e ^V(v 2 (0: t )> i anc i 
using Eq. (|114p . we find that: 



(^ 2 (0,t))-<^ 2 (0,0)) = 



d(vq) 
wo(g) 



^o(g) 

u(q) 



sin 2 w(q)t (115) 



Note, interestingly, that this result can be obtained from 
Eq. (|104p by exchanging u)o(q) and ui(q). However, when 
evaluating the integral we find that (<p 2 (0,t)) = +00, 
for all t > 0, due to the presence of infrared diver- 
gences that are not cured by the existence of a gap in 
the spectrum of Hf = H sc . Thus, we conclude that 
( e -2t4,(x^ = e -4( v 2 (o,t)) vanishes at aJi t > 0. 

The above result for the evolution of the order pa- 
rameter seems to indicate that the system apparently 
remains critical after the quench. This is conclu- 
sion is also supported by the behavior of the two- 
point correlation function of the operator e 2t ^ x,t ^: Let 
g(x,t) = (e 2 ^*)e- 2 ^(°.*)) = e « 2 c^,t) 7 w h ere C(x,t) = 
(<p(x,t)(p(0,t)) - (ip 2 (0,t)). Using Eq. © and Eq. 



uj {q) [w(q)Y 
X (1 - co$2u(q)t) (116) 



where 



C(x,0) = - ^t\(1 - cosqx)e-™ (117) 

To illustrate the above point about the apparent crit- 
icality of the non-equilibrium state, we can analyze the 
behavior of the two-point correlation function, Q[x, t), in 
two limiting cases, for t = and t — > +00. At t = 0, the 
correlation function, as obtained from Eq. (jTTTJ) , reads: 



g(x,0) = A'(a ) (2° 
V x 



(118) 



where A'(ao) depends on the short-distance cut-off ag. 
Thus, the correlations are power-law because the initial 
state is critical. In the limit where t — > +00, the part of 
the integral in Eq. (|116p containing the term cos2o;(g)t 
oscillates very rapidly and upon integration averages to 
zero. The remaining integral can be done with the help 
of tables,— yielding: 



C(x, t — > 00) —C(x, 0) 



-G 



4v 2 



1 1 

110 1/2 

(119) 

where G 2 , 2 is a Meijer function Using the asymptotic be- 
havior of the Meijer function^ we obtain: 



lim g(x,t)=B'(a )(— ) 

-^+qc \mx J 



(120) 



B'(ao) being a non-universal constant. Thus, although 
initially the system is critical and therefore correlations 
at equilibrium decay as a power law with exponent 2k 2 , 
when the system is quenched into a gapped phase (where 
equilibrium correlations exhibit an exponential decay 
characterized by a correlation length £ c w v/m), the cor- 
relations remain power-law, within the semiclassical ap- 
proximation. The exponent turns out to be smaller, equal 
to k 2 , which is half the exponent in the initial (gapless) 
state. In other words, within this approximation, it seems 
that the system keeps memory of its initial state, and 
behaves as if it was critical also after the quench. This 
behavior seems somewhat different from the results ob- 
tained for the same type of quench at the Luther-Emery 
point, where both the order parameter and the corre- 
lations for t + 00 approache a constant value, A(mao) 
(unless the non- universal amplitude A(mao) = 0, which 
seems to require some fine-tuning). Whether the differ- 
ences found here between the Luther-Emery point and 
the semi-classical approximation are due to a break-down 
of the quasi-classical approximation (which neglects the 
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existence of solitons and anti-solitons in the spectrum of 
the sine-Gordon model), or to a qualitative change in 
the dynamics as one moves away from the Luther-Emery 
point, it is not clear at the moment. To clarify this issue, 
further investigation of this issue will be required in the 
future. 



V. LONG-TIME DYNAMICS AND THE 
GENERALIZED GIBBS ENSEMBLE 

Recently, Rigol and coworkers 1 ^ observed that, at least 
for observables like the momentum distribution or the 
ground state density, the asymptotic (long-time) behav- 
ior of an integrable system following a quantum quench 
can described by adopting the maximum entropy (also 
called 'subjective') approach to Statistical Mechanics, pi- 
oneered by Jaynes i 45 ' 46 Within this approach, the equi- 
librium state of a system is described by a density ma- 
trix that extremizes the von-Neumann entropy, S — 
— Trplnp, subject to all possible constraints provided 
by the integrals of motion of the Hamiltonian of the sys- 
tem. In the case of an integrable system, if {I m } is a set 
of certain (but not all of the possible) independent inte- 
grals of motion of the system, this procedure leads to a 
'generalized' Gibbs ensemble, described by the following 
density matrix: 

PsG = J_ e -E m A m / m; (121) 

Z gG 

where Z gG — Tr A ™^ m ^phg values of the Lagrange 

multipliers, A m , must be determined from the condition 
that 

(/m)gG = Tr [p I m ] = (Im)- (122) 

where p$ describes the initial state of the system, and 
(•••) g G stands for the average taken over the general- 
ized Gibbs ensemble, Eq. (|121j) . Although pi = |$(t = 
0))(<I>(< = 0)| in the case of a pure state, as was first 
used in Ref. lid nothing prevent us from taking po to 
be an arbitrary mixed state and in particular a thermal 
state characterized by an absolute temperature, T. In 
this case, the Lagrange multipliers will depend on T or 
any other parameter that defines the initial state. 

Rigol and coworkers tested numerically the above con- 
jecture by studying the quench dynamics of a ID lattice 
gas of hard-core bosons (see Ref. [T0IIT3I for more details) . 
The question that naturally arises then is whether the 
family of integrable models studied in this work (see Eqs. 
(Trj, and their fermionic equivalences of Eq. (|rJ7|) and 
(|68[)) relax in agreement with the mentioned conjecture. 
In other words, does the average (0)(t) at long times 
relax to the value (0) gG = TrpgG O, for any of the corre- 
lation functions considered previously? In this section we 
shall illustrate some of these questions by using the Lut- 
tinger model discussed in Sect. IIIII We shall first show 
for what kind of operators the generalized Gibbs ensem- 
ble fails to reproduce their expectation values. Moreover, 



by considering the correlation function of the current op- 
erators (i.e. O = J r (x), cf. Eq IM)) , we will illustrate 
why it works. Calculations of other observables and con- 
sideration of the other models as can be found in the 
appendices. 

Let us define the generalized Gibbs ensemble for the 
Luttiger model (LM) . Since the final Hamiltonian (in the 
N = and J = sector) is diagonal in the &-boson basis, 
i.e. Hlm — J2 q ^o h v (l)\l\ b a natural choice for 

the set of integrals of motion is I m — ► I(q) — n(q) = 
W(q)b(q) for all q 7^ (a more complete version of the 
ensemble should also include N and J, but this will not 
necessary as we focus on the thermodynamic limit here) . 
Thus, for the quench from the non-interacting to the in- 
teracting state (cf. Sect. IIII B|) . where the initial state is 
|$(t = 0)) = |0), the Lagrange multipliers are determined 
by Eq. (fT2"2"]) , which yields: 

(J(g))gG = (n(q)) gG = sinh 2 P{q) = -. (123) 

Indeed, this result can be quickly established by realiz- 
ing that p gG has the same form as the density matrix of 
a peculiar canonical ensemble where the temperature on 
each eigenmode of the final Hamiltonian depends on the 
wave-vector q, that is, T(q) = liv(q)\q\/X(q). Alterna- 
tively, one can also regard it as an ensemble where the 
effective Hamiltonian that defines the Boltzmann weight 
is given by H cS /T cS = J2 q7 to H<l) n (<l)- However, it is 
worth noting that p gG is diagonal in n(q), and therefore 
it does not capture the correlations existing in the initial 
state between the q and —q modes. Mathematically, 

(n{q)n{-q)) = sinh 2 (3(q) cosh 2(5{q) ^ (n(q)n(-q)) sG 
= (n(q)) gG (n(-q)) sG = sinh 4 (5{q) (124) 

As matter of fact, since n(q)n(— q) commutes with H, we 
conclude from the above that (n(q)n(—q)} does not relax 
to the value predicted by p gG . Although this defect of p gG 
can be fixed by enlarging the set of integrals of motion to 
include the set I'(q) — n(q)n(—q) as well, we shall show 
below by explicit calculation that this is not needed for 
the calculation of the simplest correlators and observables 
in the thermodynamic limit. However, one important 
exception to this case are the squared fluctuations of the 
energy: 

<7 2 = (ffVW 2 = EMp)M9) 

x[(h(p)n(q))-(h(p))(h(q))} (125) 

which yields a 2 = 2cr 2 G = J2 q sinh2 ^l3{q)h 2 uj(q) 2 . 
Again, since the operator H 2 is conserved, a 2 violates 
the relaxation hypothesis. However, it is tempting to ar- 
gue, because that a 2 (as well as (H)) is a non-universal 
property of the LM model, this violation is less problem- 
atic than a violation in the asymptotic behavior of the 
correlation functions would be, as the latter tends to be 
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more universal. Similar results can be obtained for the 
other models considered in this work. 

Let us now consider observable correlations in the 
LM model. This requires that the operator O must 
be hermitian, which is the case for the current operator 
J r {x) — d x (f>r(x)/2iT but not for the field operator ip r (x) 
(in this case, one has to consider the momentum distri- 
bution, as we have done in Sect. lIIIBj) . In the case of the 
current operator, we shall study the following two-time 
correlation function (no time ordering is implied): 

C Jr (x,t,r) = (J r (x,t + r/2)J r (0,t-T/2)) T . (126) 

where (. . .)t stands for average other the thermal en- 
semble described by pi = e~ Hi ^ T /Z 0: with if; = H (cf. 
Eq.rjSD. Using (HU), and we obtain 



C Jr (x,t,T) = 




x {e iqx f(q, t + r/2)/*(g, t - r/2) [1 + n B {q)] 
+e* qx g*(q,t + T/2)g(q,t-T/2)n B (q) 
+e-^f* (<?, t + r/2)f(q, t - r/2)n B (q) 
+e~ iqx g(q, t + r/2)g* (q, t - r/2) [1 + n B (<?)] } , 

(127) 

being n B {q) = ( e -^o(<?)/T _ ^-1 = Vp \ q ^ the 

initial Bose distribution of modes. In the following we 
shall argue that, in the limit t — > +oo the above expres- 
sion reduces to the following correlator in the Generalized 
Gibbs ensemble: 

Cf(x,r) = Tr[ PgG (T) J r (z, r)J r (0, 0)] , (128) 

where p g o (T) is the extension to an initial thermal state 
of the Generalized Gibbs ensemble introduced above (no- 
tice that since [Hf.I(q)} = and therefore [Hi,p g Q] = 0, 
it is in principle possible to define time-dependent cor- 
relation functions on this ensemble, just as we defined 
the for standard equilibrium states ). For this (thermal) 
initial condition Q122p fixes the values of X(q) which now 
depend on f3(q) and the temperature: 

sinh 2 (3{q) [1 + n B (q)} + cosh 2 f3{q)n B {q) = A( J" • 

' (129) 

Introducing this result into the mode expansion for 
Eq. JEHU, 




x { e «K*- OT ) cosh 2 13{q) [1 + (n(q))] 

+ e lq{x+VT) sinh 2 f3{q) (n(q)) 

+ e~ iq{x - VT) cosh 2 p(q) (n(q))} 

+ e - iq(x+VT *> sinh 2 f3(q) [1 + (n(q))] , (130) 



we see that it coincides with the t — > +oo limit of (|127p . 
where the rapidly oscillating terms that depend only on 
t can be dropped 

Let us close this section with brief considera- 
tion of higher order correlation functions. In 
particular, if we consider computation objects 
like the four point current correlation function, 

{Jr{xi,ti)J r (x2,t 2 )Jr{x3,h)J r (x2,t i )) T , we WOUld 

encounter terms involving (n(q)n(—q)) , which are not 
described by the generalized Gibbs ensemble (in its 
simpler form where I(q) = n(q)). However, it is not hard 
to convince oneself that in the thermodynamic limit, 
L — > oo the contribution from such terms vanishes. This 
is because momentum conservation when computing 
( J r (x\, t\) J r (x2, t2)J r (x3, t^)J r (x2-, tifjT leaves (in most 
cases) with independent two wave numbers (out of 
four) that must be summed over. When taking the 
thermodynamic limit, the L 2 resulting from this two 
sums exactly cancels the L~ 2 factors coming from 
the mode expansions of J r (x), thus yielding a finite 
contribution. However, terms involving (n(q)n(—q)) 
require the four wave numbers to be equal, thus resulting 
in only one momentum summation, which therefore 
cannot cancel L~ 2 factor. This justifies the use a Wick's 
theorem when computing higher order correlations in 
the thermodynamic limit, but in finite size systems, 
this procedure as well as the simplest version of the 
generalized Gibbs ensemble would miss the contribution 
stemming from (n(q)n(—q)} correlations. 



VI. RELEVANCE TO EXPERIMENTS 

As we described in the introduction, cold atom sys- 
tems is the ideal arena to study quench dynamics. This 
is because they are, to a large extent, entirely isolated 
systems. Furthermore, as far as one dimensional systems 
are concerned, there are already a number of experimen- 
tal realizations, including experiments where quench dy- 
namics has been already probedi 43 ' 59 i 60 Thus, in this sec- 
tion we would like to discuss the possible experimental 
relevance of the results obtained in previous sections. As 
mentioned above, this must be done with great care, as 
our results have been obtained using effective field theory 
models that can be regarded as 'caricatures' of the Hamil- 
tonians that describe real experiments. We must empha- 
size that the situation in the case of quantum quenches, in 
particular, and of non-equilibrium dynamics, in general, 
is very different from the analysis of low-temperature 
phenomena in equilibrium. In the latter case, the ex- 
perimental relevance of effective field theories is well es- 
tablished using renormalization-group arguments and it 
has been, over the years, widely tested using a variety of 
numerical and also (when possible) analytical methods. 
By contrast, here we travel through a largely uncharted 
land, and much needs to be studied in order to achieve 
a similar level of rigor as in the equilibrium case. Thus, 
it is convenient to regard this models as 'toy models', 
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which can provide us valuable lessons and insights into 
the non-equilibrium dynamics of strongly correlated sys- 
tems. This this cautionary remarks, we can proceed to 
discuss some experimental systems to which the above 
results could of some relevance. 

The results presented in Section Unl were obtained for 
the particular case of the exactly solvable LM. The LM 
is the exactly solvable model describing the fixed point of 
a general class of interacting one-dimensional models^, 
know as Tomonaga-Luttinger liquids. This class includes 
sytems such as the one-dimensional Bose gas interact- 
ing via a Dirac-delta potential (which is solvable by the 
Bethe-ansatz^), as well as many other systems of in- 
teracting Bose gases with repulsive interactions (such as 
dipolar) or Fermi gases with both attractive and repulsive 
interactions. With the caveats of the previous section, it 
would be interesting to test the results obtained using 
the LM in one of these systems. However, the dynamics 
may be strongly modified by the fact that higher energy 
states will be also excited following a quantum quench 
such as the one described here which involves switching 
on (or off) the interactions. Such higher energy states 
are not correctly described by the LM, and the situation 
is expected to worsen in the case of short range interac- 
tions. Thus, one possible way around is to study either 
experimentally or numerically quenches where the inter- 
action does not change between too different values and 
the interactions are longer range, as the latter case pro- 
vides us with a much more faithful realization of the LM. 
One such system is a single-species dipolar ID Fermi gas 
confined to one dimension in a strongly anisotropic trap 
(e.g. — ). Since p-wave interactions are weak (away from 
a p-wave Feshbach resonace), the dominant interaction 
is dipolar, which, when the dipoles are all aligned by 
an external (electric or magnetic, depending on whether 
the dipole is electric, like in hetero-nuclear molecules, or 
magnetic, like in Chromium). When confined to ID, the 
dipolar interaction between the atoms can be approxi- 
mated by the potential: 

^• 9) = ipTlS^ (131) 

where D is the dipolar momentum of the atoms, 9 is the 
angle subtended by the direction of the atomic motion 
and the polarizing field, and X(6) = (1 — 3cos#). Since 
in this case 32(9) = 54(9) oc X(9), a sudden change in 
the interactions can be produced by a sudden change in 
alignment of the field with the direction of motion, that 
is, a change in 9. In particular, a change in 9 away from 
the value 9 m = cos~ 1 (|) would produce lead to a sud- 
den switching of the interactions amongst the Fermions. 9 
At zero temperature, the momentum distribution /(p, t) 
(which can be probed by time of flight measurements) 
following the quench into the interacting system would 

evolve as described in Sect. IIII Bl (cf. Fig. [1]), with the 

2 

discontinuity at the Fermi level dying out as t 1 . How- 
ever, currently atomic gases are produced at tempera- 
tures T ~ 10% to 20% of the Fermi energy, and this 



would complicate the observation of this effect. If much 
lower temperatures could be reached in experiments, so 
that the application of effective field theory is much more 
reliable, we expect that in a time of the order of H/T the 
quenched dipolar gas reaches a stationary state charac- 
terized by a momentum distribution that differs from the 
thermal one. However, the calculations of f(p,t) pre- 
sented in Sect. [HTbI (cf. Figs. H and |3J) show that the 
differences between the non-equilibrium and equilibrium 
results in the stationary state may be well below the cur- 
rent experimental resolution. Alternatively, instead of 
measuring the momentum distribution, one can try to de- 
termine the non-equilibrium exponents measuring noise 
correlations in the time-of-flight images 61 or through in- 
terferometry.— 

Finally, we shall remark that the sine-Gordon model is 
an effective field-theory description of the Mott insulator 
to superfluid transition (MI to SF) in 1D 55 i 56 or, when 
the field <p is interpreted as the (relative) phase of two 
coupled ID Bose gas, it describes the Josephson coupling 
of two ID Bose gasesi 39 i 56 Thus, the results presented 
here could be of some relevance when interpreting the 
evolution of such systems following a quantum quench. 
In particular, in the case of the MI to SF transition, 
the order parameter e~ 2i ^ x ' describes the time evolution 
oscillatory part of the density in the system, whereas the 
correlations of the order parameter are related to the time 
evolution of the static structure factor. 
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APPENDIX A: DETAILS OF THE 
CALCULATION OF THE ONE-PARTICLE 
DENSITY MATRIX IN THE LUTTINGER 
MODEL 

In this Appendix, we shall provide the details of the 
calculation of non-equilibrium one-particle density ma- 
trix: 

C^(x,t) = (e iH ^ h 4(x)i> r (0)e- mt/h ), (Al) 
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To this end, the formula (p8|) is used. In normal ordered 
form: 



tp a (x) = 



TX / L 



(A2) 



where the normal order is defined: 

. s a i(p a (x) . pitpct pS^irixNct s a i^ a (x) s a i<& a (x) 

• c • — O O o o ■ 

The boson field $q,(x) is given by Eq. ([30| . Hence, 

: e-^ix) ■ ■ e +i4>A0) . _ e -2wixN r /L e [<£r(x),$U0)} 



(A3) 



where we have used the identity e e 



,[A,B] p B p A 



which holds provided [A, B] is a c-number. Using that 
(ao — * + is the short-distance cut-off): 



2tt 



q>0 

= - In [l - e -^ /L e 2i^ x /L : 

we arrive at the following expression for C 1 f, r (x 1 t): 
C^{x,t) = G^(x) ^-iFH^) e -iFAx,t) }j 

where 

G^(x) 



1 



2L sin [f (x + m)] 
F r (as,t) = e ifft/ri [* r (x) - $ r (0)] e^*^ 

<3>0 Vy 7 



(A5) 
(A6) 

(A7) 

(A8) 
(A9) 
(A10) 



X [f(q,t)b(q)+g*(q,t)tf(-q)] . 

(All) 

To derive the last expression we have used Eq. |7]). Em- 
ploying the identities e A e B — e^ A ' B ^ 2 e A+B (provided 
[A,B] is a c-number) and that (e D ) = e (D)+±((D-(D)) 2 ) | 
we obtain: 

O/vM - GjpJ^e-^Cx.tJM*,*))^ (A12) 

where we have used that ([F r +(x, <)] 2 ) = (F 2 {x,t)) = 
because ([& f (g)] 2 ) = (Hq)} 2 ) = (tf(-q)b(q)) = 0, and 
since the commutator [F$ (x, £), F r (x, t)] is a c-number, it 
can be safely replaced by ([F^(x, t), F r (x, t))). Note that 
for t — F r (x,t — 0) contains only b(q) and thus the 
average (F}(x, i)F r (x, t)) = Q at T = 0. In Eq. (|AT2)) the 
exponent can be expanded to yield: 



q>0 



(F}(x,t)F r (x,t))=J2(^)e 



qL 



x [|/(9, t)\ 2 n (q) + \g(q, t)\ 2 (n (q) + 1)] , (A13) 

being no((j) = (e^' 9 ' — l) -1 the distribution of Tomonaga 
bosons in the initial state (which has been assumed to be 
a mixed thermal state), and A = hvp/T is the thermal 
correlation length. We next evaluate explicitly the above 
result in several limiting cases. 



1. Zero temperature and finite length 

Let us now consider the T = limit, where no(q) = 0, 
and thus, using © and ©, Eq. (|A13[) simplifies to: 

(FHx,t)F r (x,t)) T=0 = £ (^L) e-« a ° sinh 2 [2/?(g)] 



q>0 



x (1 - cosgx) [1 - cos2v(q)\q\t}. (A14) 

To make further progress, we shall assume that 
sinh 2/3(q) = ^e^^ qRa ^ 2 where Rq is the range of the inter- 
action. Furthermore, we shall replace v(q) by v = i>(0) 70 , 
what allows us to safely take the limit ao — > 0+. Next, 
in order to simplify the computation, we introduce the 
quantity 



(A15) 



which can be readily computed to give 



£ r {z) = - In -d(z + iRo\L) 

- 1j 



L 



-In 2, (A16) 



where d(z\L) — L \ sin(7rz/L)|/7r is the cord function. Us- 
ing this result into Eq. (|A14j) . yields the following ex- 
pression for the one-particle density matrix: 



C 4 , r {x,t>Q\L) = Gf\x\L) 



d{iR \L) 



d{x + iR \L) 



d(x + 2vt + iR Q \L)d(x ~ 2vt + iR \L) 



d(2vt + iR \L)d(-2vt + iR a \L) 



7 2 /2 



(A17) 



Taking into account that Rq/L <C 1 and, in the scaling 
limit, one recovers the result quoted in the main text, 
Eq. (32]). 



2. Thermodynamic limit and finite temperature 



We next consider Eq. (|A12[ ) for L — ► oo and finite tem- 
perature, T. Equation (|A13|) can be recast as: 

{F$(x,t)F r (x,t)) T = (F r (x,t)F r (x,t)) T=0 

+ H(x) + g(x,t) (A18) 

where we have introduced the following functions: 

H(x)=2 [ — e- qa {l- cos qx)n (q), (A19) 
Jo 1 

g(x, t) = 2 7 2 / ^- e ~ qRo (1 - cos qx) (A20) 

Jo q 

[l-cos(2vqt)]n (q), (A21) 

which hold in the thermodynamic limit and upon replac- 
ing v(q) by v = v(q = 0) and sinh2/3(g) = ^e~^ qRa ^ 2 
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as we did in the previous section. We next define the 
function 



which can be evaluated to yield: 57 

r(l + A"V) 



j(u; r) = 2 In 



r[l + A-!(r 



(A22) 



(A23) 



where T(z) is the Gamma function. In the limit where 
r <C u, and using that 5 ^ T(z)T(l — z) = it/ sin(7rz), the 
above expression reduces to 



g{u;r) = - In 



dh{ir\T) 



-In 



u + ir 



dh(u + ir\T) 
In the previous expression we have defined: 

dh(z\T) = - Isinh (7rA _1 z) I . 



(A24) 



(A25) 



Combining this result with Eqs. (|A18[) - (|A21|) and (IA12I) . 
it is seen that the second term in Eq. (IA24j) ex- 
actly cancels the contributions from G^°'(x\L) and 
{F r (x,t)F r (x,t))T=o m the thermodynamic limit, and 
therefore, 



C^ r (x,t\T) = G^(x\T) 



dh(iR \T) 



dh(x + iR \T) 
dh{x + 2vt + iRo\T)dh(x - 2vt + iR^Ty 1 ' '' /2 



dh{2vt + iR \T)dh(-2vt + iR \T) 



where 



Gf)(z|T) 



ttA" 



27r sinh [tt\ 1 (x + iao) 



(A26) 



(A27) 



We note that the result of Eq. (|A26p can be obtained from 
Eq. (|A17[) upon making the replacement L sui(irL~ 1 x) jit 
by Asinh(7rA a;). Taking into account that Rq/L, and 
in the scaling limit, we retrieve the result quoted in the 
main text, Eq. 



Using the bosonization identity, Eq. (|28|) . we can write 
the expression as follows: 



C f r {x) = Gf\x) (: e -*[M-)-M0)] :)gG 
= G(°)(i) (e-^^e-^.)^ 
Taking into account that 



F r (x) = $ r (x) - $ r (0) 



E 

9>0 



2tt\ 



1/2 



(B2) 
(B3) 

(B4) 
(B5) 



[cosh (3{q)a{q) - sinh/3(q)a t (-( ? )] . (B6) 



The expression for (x) can be easily computed by 
using the trick of regarding p g c as a canonical ensemble 
with g-dependent temperature. Thus, following the same 
steps as in the previous section we arrive at: 



Cf r {x)=G?\x) e-<ftWM*))&. 

Given that 

(F}(x)F r (x)) gG = sinh 2 2f3 [V r (0) - V r (x)} 
where 



V r (x)=Re^(^z) 6 " 



qRo e iqx 



= — In sin —(x + iRq) 
L 



-ln2- 



irRn 



Hence, taking the thermodynamic limit 



cT{x) 



2ix(x + ia) 



Rc. 



X 



(B7) 
(B8) 

(B9) 
(BIO) 

(Bll) 



Thus we see that one recovers the same results as 
lim t ^ +00 Gp T (x,t), Eq. 



APPENDIX C: THE SINE-GORDON MODEL 
AND THE GENERALIZED GIBBS MODEL 

1. Quench from the gapped to the phase at the 
Luther-Emery point 



APPENDIX B: ONE-BODY DENSITY MATRIX 
OF THE LUTTINGER MODEL IN THE 
GENERALIZED GIBBS ENSEMBLE 

Next we take up the calculation of the one-body den- 
sity matrix in the generalized Gibbs ensemble for the 
Luttinger model discussed in Sect. El That is, we shall 
evaluate the expression at T = 0. 

Cf r (x) = Tr[p sG 4(x)M0)] (Bl) 



In this case, the type evolution of the system is per- 
formed by Hq (cf. Eq. I66[) . which is diagonal in the op- 
erators n a (p) =: ipa(p)4>a(p) ■ (a = r, 1). Thus, the 
generalized Gibbs ensemble can defined by the follow- 
ing set of integrals of motion /,„ — * I a (p) = n a (p). We 
see immediately that the fact this ensemble is diagonal 
in n^{p) means that the order parameter, {e~ 2ltp ^) g Q — 
{tp r {x)ipi(x)) g G = 0, which agrees with the t — ► +oo limit 
of the order parameter, which was shown in Sect llV^l to 
exhibit an exponential decay to zero. However, the two- 
point correlator of e 2lv ^ has a non-vanishing limit for 
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t — > +00. Thus, our main concern here will be the calcu- 
lation of the correlation function: 



(e -2 M x) e 2M0) )gG = ($(sM(*)$(0ty r (0)> 



gG 



E 



Pl,P2,P3,P4 



a i(pi-p 2 )x 
Z2 



(lpl(j3l)lpl (P2 ) ^/ (P3 ) (P4 ) ) gG 



(CI) 



Since the ensemble is diagonal in the chirality index, a, 
as well as momentum, p, we evaluation of the above ex- 
pression can be carried out by noting that: 



Tr e~ Cp)'.' (Ptyt (p)^ Q (p) 



e A«(p) +1' 



(C2) 



where the Lagrange multipliers A(g) can be related to the 
values of the same expectation values in the initial states 
by imposing their conservation, that is, 



(V4(p)Vv(p)) g G 



e M P ) + l (C3) 
($(*#l(p)) =cos 2 0(p) (C4) 
1 



eM« + 1 ( C5 ) 
(^(p)^(p)>=sin 2 0(p). (C6) 



Hence, 



(V'J(Pl)V'i(P2)'0i t (P3)V'r(P4))gG = (^(pi)V'r(P4))gG 

= 5 Pl , Pi 5p 2 , P3 sin 2 0(pi) (1 - cos 2 0(p 2 )) 

= 8 Pu pJp2,P3 sin2 #(Pi) sin2 ^(Pa). (C7) 
Introducing the last expression into Eq. (|C1[) yields: 



(e -2 lV (x) e 2 4V (0) )gG 



- 5^ e apa; sin 2 6{p) 



(C8) 



and using that sin 2 6(p) = (1 — cos 2 2(9(p))/2 and 
cos 29 (p) = ojo (p) I \Ju>q (p) + to 2 , we find (for x 7^ 0), 



-2iy( 



^ (0) >gG = (^Y 

V 2tTV I 



, (C9) 



which equals to the t — > +00 limit of Eq. ([85 



2. Quench from the gapless to the gapped phase at 
the Luther-Emery point 

In this case the initial state is the gapless ground state 
of Hq, Eq. 1661 whereas the Hamiltonian that performs 



the time evolution has a gap in the spectrum and it is 
diagonal in the basis of the ipv(p) and ip c {p) Fermi op- 
erators (cf. Eq. [71]). Therefore, the conserved quantities 
are 

I v {p)=n v (p)=ipl{p)^{p), (CIO) 
Ic(p) = n c (p)^ c (p)Mp) (Cll) 

The associated Lagrange multipliers (at zero tempera- 
ture), X v (p) and \ c (p) can be obtained upon equating 
(IvAp)) S g = (¥(0)|I OlC (p)|*(0)>. This yields: 



(IviP)), 



(Ic(p)), 



1 



- - (C12) 

e Mp) + l v ' 

= §{-p) sin 2 9(p) + 'dip) cos 2 9{p) (C13) 

= -T-ri (C14) 

e K(p) + 1 y J 

= d(-p) cos 2 Oip) + dip) sin 2 0(p), (C15) 



where flip) denotes the step function. Using these ex- 
pressions we next proceed to compute the expectation 
values of the following observables: 



a. Order parameter 
We start by computing the order parameter, 
(e- 2i ^) gG = (tfnWM*)) (C16) 

= iE Sin20 (P) [< J c(P)) g G - (Ivi P )) gG ] , (C17) 

p 

and upon using Eqs. (|C13p and (|C15p . 



-^sin20(p) cos 29 (p) 



p>0 



-pao 



(C18) 

(C19) 
(C20) 



dp mujo(p) 
1 2ir uj 2 (p) 
= A(mao) 

where we have used that cos 29- p = - 
is the non- universal constant introduced in Sect. IIVB 21 
This result agrees with the one obtained in Sect. IIV B 21 
for the order parameter in the limit t — > +00. 



cos 29 p \ A(mao) 



b. Two-point correlation function 

We next consider the two-point correlator of the order 
parameter, namely 



(e 2M,) e 2M0) )gG = j_ J- 



i{pi-pi)x 1 i\ 



Pl,P2,P3,P4. 

x V , /(P3)V , r(>4))gG- 



{lpl{pi)tpl{p 2 ) 



(C21) 
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The calculation of the average in this case is a bit more 
involved, but it can be performed by resorting to a factor- 
ization akin to Wick's theorem. This is applicable only 
in the thermodynamic limit, as it neglects terms the four 
momenta of the above expectation value coincide. These 
terms yield contributions of 0(1/ L) compared the oth- 
ers. When factorizing as dictated by Wick's theorem, the 
only non- vanishing terms are: 



(V4 (Pi ) i>l (P2 m(P3)lpr(Pi)) gG 
X (ipt (Pi ) Ipr (P4 ) ) gG (ip j (P3 ) i>l (P2 ) ) gG 



(C22) 



&P1P2&P4.P3 (4(Pi)MP2)) gG (ri(P3)MP±)) s G- (C23) 



Upon using 

(4(p)Mp)) f 



-d(p) sin 2 20 (p) 



+ 0(-p)(l- ^ sm 2 26(p) 



(V, + &#j(p)>gG = ^(-P) sin 2 2%) 



■0(p) ( l--sin 2 2%) 



(C24) 



(C25) 



(il)\(p)ip r (p)) &G = (ipi(p)^i(p))gQ 



= --sin2(9(p)cos26>(p)sgn(p), (C26) 

the average over the generalized Gibbs ensemble of the 
four Fermi fields on the right hand-side of Eq. (|C21[) can 
be computed and yields the following expression for the 
two-point correlation function (up to terms of 0(1/L 2 )): 



and coincides with the t — > +oo limit of the second term 
in the right hand-side of Eg. [Ml in Sect. IIV A 21 (the func- 
tion H(x,t) is defined in Eq. [ 



3. Quench from the gapless to the gapped phase in 
the semi-classical approximation 

In this case Hamiltonian performing the time evolu- 
tion is gapless, and therefore, diagonal in the 6-operators. 
Hence, the conserved quantities are 



I(q) = tf(q)b(q) 



(C29) 



The Lagrange multipliers of the corresponding general- 
ized Gibbs density matrix are fixed from the condition: 



1 



e A(<r) - i 
= sinh 2 f3(q), 



($(0)|6t(«)6(g)|$(0)) 



(C30) 
(C31) 



where (3(q) is defined by Eq. (llOlj) Hence, using this 
result we next proceed to compute the order param- 
eter and the two-point correlation function. We first 
note that the order parameter vanishes in the general- 
ized Gibbs ensemble since (e- 2 ^ x )) gG = e ~ 2l ^°^ G 

and (0 2 (O)) g G = x(^ 2 (0)) is divergent in the L — > +oo 
limit (see below). This agrees with the result found in 
Sect. IIVB 11 where it was found that the order parame- 
ter decays exponentially in time. Thus, in what follows 
we shall be concerned with the the two-point correlation 
function. 



(e 2 ^V 2i ^°)) gG 



ip.r 



§(j>) sin 2 26 (p) 



+ 0(-p) ( 1 - - sin" 26 p 



-- Y]sin20(p) cos 29 (p) 

p>0 



(C27) 

The last term in the above expression in just 
(e 2tv{x) ) gG (e~ 2lv{0) ) g G (cf. Eq. lCi"8|) . whereas the first 
term in the right hand-side can be written as 



(-P) + 2 s S n (p) sin2 2e v 



p>0 



-ipx 



- ^2 sin P x 



p>0 



UJ 2 (p) 



in 

p>0 



lim H.(x,t) 

t — > + oo 



(C28) 



a. Two-point correlation function 



Since (e^^e 2 ^ ))^ = e -4 CEG (*), where 
C° G (x) = (<p(x)(p(p))gB - {ip 2 (0)) gG . In order to ob- 
tain this correlator, we introduce the Fourier expansion 
of ip(x) (neglecting the zero-mode contribution), 



2nv 



2 ^ \u (q)L 



1/2 



e^[b(q) + b\-q)\, (C32) 



into the expectation value, and using (|C31[) to evaluate 
the averages in the generalized Gibbs ensemble, we find 
that, in the thermodynamic limit, 



(ip(x)(p(0)) gG 



d(vq) 



cosqx cosh 2 (3(q), (C33) 



l 4w (g) 

and therefore, 

V G (x) = M*M0)) gG - <^ 2 (0)) gG , (C34) 



C 9G (x) = 



d(vq) 
o vo(q) 

2 /> + oo 



cosh 2/3(q) (1 - cos qx) (C35) 



C(x,0)-— I d{ - Vq) J l- cosqx) (C36) 

4 J w (<z)M<7)] 2 
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where C(x, 0)=C(x,t = 0) is defined in Eq. (fTTOf . Upon 
comparing the last result with Eq. p09p in the limit 
where t — > +00, we see they are identical. 



4. Quench from the gapless to a gapped phase 



e 2 



<v 2 (o)) gGj an d (^(O))^ is logarithmically di- 
vergent in the thermodynamic limit (see expressions 
below), the find that ( e - lK ^) gG = 0. This r esult is 
in agreement with the one found in Sect. IIV B 21 for the 
order parameter. 



In this case the Hamiltonian that performs the time 
evolution is gapped, whereas the initial state is gapless. 
Thus, differently from the previous case, the Hamilto- 
nian that performs the evolution is diagonal in the a- 
operators, and therefore, the conserved quantities are 
I(q) = a^(q)a(q). The corresponding Lagrange (at zero 
temperature) are fixed from the condition: 

™>sG = J = (mW(q)a(q)\m) (C37) 



sinh 2 (3(q), 



(C38) 



where (3(q) is given by Eq. ()101jl . 

In order to obtain the one and two-point correlation 
functions of e 2 ^ 2 ^ = e 2lKip ^ we first need to write the 
tp(x) field in terms of the a-operators. Upon using the 
canonical transformation Eq. (|3j): 



ip(x) 



2?r 



2 V \u(q)L 



1/2 



[0(9)+ at (C39) 



Hence, since (e 



)gG 



-ZKtp 



(*)) 



a. Two-point correlation function 

Next we consider the two-point correlation function 
of the same operator, namely (e -2 *^ 3 ^ 2 "^ ))^ = 

e-^ csC H where C^z) = (^M0)} gG - (^ 2 (0)) gG . 
We first obtain: 



MxMo)> gG = 



djvq) 
M<z) 



cos g£ cosh 2/3(g). (C40) 



Hence, 
C gG (a;) 



djvq) 
M<?) 



cosh /?(<?) (1 - cos qx) (C41) 



C(x,0) 



2 /-+oo 



d(vq) 



4 7 ^o(g)M<?)] : 



r (l — cos qx) 
(C42) 



where C(x,0) is defined in Eq. (I117p . The latter result 
agrees with Eq. ()116|) in the t — > +00 limit. 
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Implicit in this discussion it is the fact we are attempting 
at a description of the system as a whole, not separating 
its degrees of freedom into a 'system' and a 'reservoir'. We 
believe that this point of view may be the more appropriate 
when discussing cold atomic systems, given that they are 
very weakly coupled to the environment. 
This question is also related to the problems concerning 
the applicability of the maximum entropy approac h 45 ' 46 to 
Statistical Mechanics. See for instance the critique by S. K. 
Ma in Ref. l63l.. However, the maximum entropy approach 
is advocated by R. Balliant.— 

In a Luttinger liquid, the correlator describes the fluc- 
tuations of wave number close to 2mltF (where &f is the 
Fermi momentum) of the density-correlation functio n 54,56 . 
In the presence of a periodic potential of periodicity equal 
to 2mkp the system may become an insulator.— The 
power-law behavior exhibited at zero temperature, and in 
the thermodynamic limit) is a consequence of the gapless 
spectrum and the absence of long range order in the den- 
sity. In the insulating (i.e. gapped) phase, this correlation 
function decays to a non-zero constant at long distances, 
which is a consequence of the existence of long ranger order 
in the density at wave number 2mfci?. Similarly, in equi- 
librium Cg measures the phase fluctuations, and exhibits 
a power-law, reflecting the absence of long range order in 
the phase. However, in the Luttinger model of interest for 
us here, terms with m > are absent from the density 
operator, which is given by p(x) — J r (x) + J; (a:).— 
We shall encounter a similar situation in Sect. IIVI when 
analyzing quenches at T = from a non-critical (that is, 
gapped) into a critical (that is, gapless) state. In that case, 
the role of A will be played by the correlation length of the 
system that is determined by the (inverse of the) energy 
gap in the initial state. 

This case was studied earlier by Calabrese and Cardy in 
Ref. HH although not as limit of the sine-Gordon model. 
We assume implicitly that both f3(q) and v(q) are not sin- 
gular at q = 0. 



